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Abstract 



Within the framework of unified gauge models, interactions responsible 
for neutrino masses can also provide mechanisms for nucleon instability. We 
discuss their implications concretely in the light of recent results on neutrino 
oscillation from the SuperKamiokande collaboration. We construct a pre- 
dictive S*O(10)-based framework that describes the masses and mixing of all 
quarks and leptons. An overconstrained global fit is obtained, that makes five 
successful predictions for quarks and charged leptons. The same description 
provides agreement with the SuperK results on atmospheric neutrinos and 
supports a small-angle MSW mechanism. We find that current limits on nu- 
cleon stability put significant stress on the framework. Further, a distinctive 
feature of the S'O(IO) model developed here is the likely prominence of the 
fi + K° mode in addition to the VK + mode of proton decay. Thus improved 
searches in these channels for proton decay will either turn up events, or force 
us outside this circle of ideas. 

I. INTRODUCTION 

Recent SuperKamiokande observations on atmospheric neutrinos QXJ] establish the oscil- 
lation of to v T with a mass splitting 5m 2 ~ (10~ 2 to 10~ 3 ) eV 2 and an oscillation angle 
sin 2 26*°^ = (0.82 — 1.0). To be more precise, the observations do not directly exclude oscil- 
lation into some other ux, as long as it is not to z/ e , but Occam's razor and the framework 
adopted in this paper suggest X = r, and we shall assume that in what follows, without 
further comment. These observations clearly require new physics beyond what is usually 
contemplated in the Standard Model. 

As we shall presently discuss, a tau neutrino mass consistent with the observed oscilla- 
tions fits extremely naturally into the framework supplied by unified gauge theories of the 



strong, weak, and electromagnetic interactions that includes the symmetry structure 
G 224 = SU(2) L x SU{2) R x SU(4) C @, the minimal such symmetry being 50(10) f|. In this 
framework, the low-energy degrees of freedom need be no larger than the Standard Model. 
Neutrino masses appear as effective nonrenormalizable (dimension 5) operators. 

Unified gauge theories were already very impressive on other grounds. They combine 
the scattered multiplets of the Standard Model (five per family) into a significantly smaller 
number (two for SU(5), one for 50(10)). They rationalize the otherwise bizarre-looking 
hypercharge assignments in the Standard Model |2||| . Finally, especially in their supersym- 
metric version ||, they account quantitatively for the relative values of the strong, weak, 
and electromagnetic couplings 

This last feat is accomplished by renormalizing the separate couplings down from a single 
common value at a unification scale, taking into account the effects of vacuum polarization 
due to virtual particles, down to the much lower mass scales at which they are observed 
experimentally. A by-product of this overconstrained, and singularly successful, calculation, 
is to identify the mass scale at which the unified symmetry is broken, to be M\j ~ 2 x 10 16 
GeVe. 

This value is interesting in several respects. First, from data and concepts purely internal 
to gauge theories of particle interactions, it brings us to the threshold of the fundamental 
scale of quantum gravity, namely the Planck mass 2 x 10 18 GeV (in rational units). Reading 
it the other way, by demanding unification, allowing for both the classical power-law running 
of the gravitational coupling and the quantum logarithmic running of gauge couplings, we 
obtain a roughly accurate calculation of the observed strength of gravity. 

Second, it sets the scale for phenomena directly associated with unification but for- 
bidden in the Standard Model, notably nucleon decay and neutrino masses. Prior to the 
Super Kamiokande observations, the main phenomenological virtues of the large value of the 
unification mass scale were its negative implications. It explained why nucleon decay is 
rare, and neutrino masses are small, although both are almost inevitable consequences of 
unification. Now the scale can also be positively identified, at least semi-quantitatively. 

Indeed, any unification based on G 224 @ requires the existence of right-handed neutrinos 
v R . When O224 is embedded in 50(10), v R fills out, together with the 15 left-handed quark 
and lepton fields in each Standard Model family, the 16 dimensional spinor representation 
of 50(10). The v R are Standard Model singlets, so that they can, and generically will, 
acquire large Majorana masses at the scale where unified 50(10) symmetry breaks to the 
Standard Model SU(3) x SU{2) x £7(1). The ordinary, left-handed neutrinos couple to these 
u R much as ordinary quarks and leptons couple to their right-handed partners, through 
SU(2) x £7(1) non-singlet Higgs fields. For the quarks and leptons, condensation of those 
Higgs fields transforms such interactions directly into mass terms. For neutrinos the effect 
of this condensation is slightly more involved. As mentioned, the v R have an independent, 
and much larger, source of mass. As a result, through the "see-saw" mechanism |J, the 
effective masses for the left-handed neutrinos, acquired through their virtual transitions into 
v R and back, are predictably tiny. 

In this paper we do two things. First, we flesh out with quantitative detail the rough 
picture just sketched. Its most straightforward embodiment leads to the hierarchical pattern 
m Ue « m Ufl « m UT , and to a value of m„ T very consistent with the SuperK observations, 
interpreted as — u T oscillations. We demonstrate that the large mixing angle observed 
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does not force us to swerve from this straightforward direction. Indeed, it can arise rather 
plausibly in the context of ideas which have been applied successfully to understanding 
quark masses and mixings. Motivated by the success of the circle of ideas mentioned above, 
we insist that the pattern of neutrino masses and mixings should be discussed together 
with those of the quarks and charged leptons, and not in isolation. A simple and predictive 
50 (lO)-based mass structure that describes the observed masses and mixings of all fermions 
including those of the neutrinos will be presented and analyzed. We thus demonstrate by 
example how the large — v T oscillation angle can be obtained quite naturally along with a 
large hierarchy in the — v T masses. (This is in contrast to several recent attempts where 
such a large oscillation angle is explained as a consequence of the near degeneracy of the 
v^ — Vr system. The bulk of the papers in Ref. P,|lO| belong to this category.) In this scheme, 
v e — oscillation drives the small angle MSW explanation of the solar neutrino puzzle. 

Second, we revisit a previously noted link between neutrino masses and nucleon decay, in 
the framework of supersymmetric SO (10) unified models ||12|| . Previously, motivated in part 
by possible cosmo logical indications for a hot dark matter component, we used numerical 
estimates for m Vr ~ 1 eV, considerably larger than are now favored by the SuperK result 
(~ 1/20 eV). This amounts to an increase in the Majorana mass of compared to previous 
work, and correspondingly an increase in the strength of the neutrino mass related d — 5 
proton decay rate. Another important change is caused by the large — v T oscillation 
angle suggested by the SuperK result. With hierarchical neutrino masses, we will argue, 
their result strongly suggests substantial mixing in the charged lepton (\i — r) sector. That, 
in turn, affects the strength and the flavor structure not only of the neutrino related, but 
also of the standard d = 5 proton decay operators induced by the exchange of color triplet 
partners of the electroweak Higgs doublets These adjustments in our expectations for 
proton decay turn out to be quite significant quantitatively. They considerably heighten the 
tension around nucleon decay: either it is accessible, or the framework fails. 

As an accompanying major result, we observe that, in contrast to the minimal supersym- 
metric SU(5) model for which the charged lepton mode p — > fi + K° has a negligible branching 
ratio (~ 10~ 3 , see text), for the 5*0(10) model developed here the corresponding branching 
ratio should typically lie in the range of 20 to even 50%, if the relevant hadronic matrix el- 
ements have comparable magnitudes. This becomes possible because of the presence of the 
new d = 5 proton decay operators mentioned above, which are related to neutrino masses. 
Thus the fi + K° mode, if observed, would provide an indication in favor of an 5*0(10) model 
of the masses and mixings of all fermions including neutrinos, as developed here. 

This paper is organized as follows. In Sec. II, we discuss the scale of new physics implied 
by the SuperK observations. In Sec. Ill we describe a caricature model that accommodates 
large neutrino oscillation angle as suggested by SuperK without assuming neutrino mass 
degeneracy. Sec. IV is devoted to a more ambitious 5O(10) model that accounts for the 
masses of second and third generation quarks and leptons including the large neutrino oscil- 
lation angle. In Sec. V we suggest, by way of example, a predictive way to incorporate the 
first family fermions into the 5O(10) scheme that retains the success of Sec. IV, leading to a 
total of eight successful predictions for the masses and the mixings of the fermions including 
the neutrinos, and supports a small angle MSW mechanism. In Sec. VI we discuss the 
issue of proton decay in the context of neutrino masses. Four Appendices (A,B,C and D) 
contain relevant technical details of our proton decay calculations including unification scale 
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threshold corrections to a^imz)- Finally, a summary of our results and some concluding 
remarks are given in Sec. VIlQ. 



II. m Vt and the unification scale 

Using the degrees of freedom of the Standard Model, small Majorana masses for neutrinos 
arise from dimension-5 operators in the form [fTH 



= M 

where Lj = (h / i ,£ i ) T denote the lepton doublets and = (0 + ,0°) T the Higgs doublet of the 
Standard Model. Interpreting the SuperK result as a measure of m UT = (1/30 — 1/10) eV, 
momentarily ignoring mixing, and using (0°) = 246 GeV, we find M/A33 = (6 — 18) x 10 14 
GeV. 

According to the seesaw mechanism, and again putting off the question of mixing, the 
tau neutrino mass is given as 

^ = "-iff W 

where m D is the Dirac mass that the neutrino would acquire in the absence of the large 
Majorana mass of the right-handed neutrino, and M R is the value of this Majorana mass. 

If one assumes, within S'O(IO), that the Dirac masses of the third family are dominated 
by a contribution from a fundamental Higgs 10h condensate then one obtains the relation 
m T (Mu) = mfe(Mu) between the masses of the tau lepton and bottom quark at the unifica- 
tion scale, which is known to be successful This suggests that the third family fermions 



get their masses primarily from a 10h condensate through a Yukawa coupling 16316310k- 
(16j, i = 1 — 3 denotes the three generations of fermions.) This hypothesis entails the 
relation 



m^(Mu) « m t (M v ) » (100 - 120) GeV . (3) 

(This numerical estimate is valid for most values of the MSSM parameters, so long as tan j3 
is neither too large (> 30) nor too small (< 2).) Combining this with the seesaw formula 
and the SuperK measurement^, one obtains M R « (1 — 3) x 10 14 GeV. 



Preliminary results of this investigation were announced at summer conferences [14|. 

2 As is well known, for very low values of tan/3 < 2, which corresponds to a large v T Yukawa 
coupling, renormalization group effects lead to significant increase in m VT , by more than 40%, 
as the running momentum (fj,) is lowered from M\j to 1 GeV [17|. Such effects however become 
progressively smaller (< 12%) for tan/3 > 3. Barring very low values of tan/3 < 2 which seem to 
be disfavored (see remarks later), such renormalization effects on m Vr are thus small. Henceforth 
we neglect them. 
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In SO (10), the Majorana masses of right-handed neutrinos can be generated using either 
the five-index self-dual antisymmetric tensor 126h, or a bilinear product of the spinorial 
Higgs 16h- The relevant interactions are the renormalizable interaction /jjl6il6jl26H or 
the effective nonrenormalizable interaction fijl6il6jl6ul^ii/M, respectively. In terms of 
the SU(2)l x SU(2)r x SU(4)c subgroup the relevant multiplets transform as (1,3,10) h or 
(l,2,4)/f respectively 

If the 126h is used to induce the Majorana mass of v^, then with (126h) ~ M\j we 
require the Yukawa coupling / 33 ~ 10~ 2 . Considering that it was constructed by balancing 
quantities of vastly different magnitudes, the nearness of this coupling to the 'natural' value 
unity is encouraging. 

Still more interesting is the situation that arises if we employ the 16h- In this case, as 
mentioned just above, we require an effective nonrenormalizable interaction. Such an inter- 
action could well arise through the exchange of superheavy states associated with quantum 
gravity. Then using (l6 H ) ~ M v ss 2 x 10 16 GeV and M « M Planck ss 2 x 10 18 GeV we find 

« fJ^^ « (/ 33 ) x 2 x 10 14 GeV . (4) 

This is nicely consistent with the required value of M^, for / 33 close to unity !0 

Many of the considerations that follow do not depend on which of these alternatives, 
16h or 126#, is chosen. Motivated partly by the foregoing numerology, and partly by some 
suggestions from higher symmetry schemes and string theory ||18|| , we will mainly discuss 
models in which the pair (16h, 16h) is used to break (B — L). Let us note, however, that 
the 126# does have some advantages. Specifically, its couplings are renormalizable, 
and its vacuum expectation value violates (B — L) by two units, and thus conserves an 
R parity automatically |jn|. The latter property is important for eliminating catastrophic 
- dimension 4 - sources of proton decay. With the SO (10) spinor condensates, we must 
postulate a suitable Z 2 symmetry for this purpose separately. 

III. LARGE fa - v T ) OSCILLATION ANGLE WITH HIERARCHICAL MASSES 

Based on its measurements of atmospheric cosmic ray neutrino oscillations, the SuperK 
group estimates a large oscillation angle sin 2 2#°^ c = (0.82 — 1) If we compare this to the 
analogous angle for quarks, sin 2 2# C 6 « 4|V c f,| 2 ^ 6 x 1(V 3 , a challenge arises. How are we 
to understand the enormous difference between these two mixings, in a framework where 
quarks and leptons are unified? 

One widely considered possibility is to propose that the two relevant neutrino flavor 
eigenstates are nearly degenerate PJItj. Then a small perturbation that lifts this degeneracy 
will induce maximal mixing. Such behavior is familiar from the K° — K° system. But while 
there is a well-established fundamental symmetry (CPT) which guarantees exact degeneracy 



3 The effects of neutrino mixing and of possible choice of M = M s t r i ng ~ 4 x 10 17 GeV (instead of 
M = Mpi anc k) on are discussed later. 
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of K° and K°, no such symmetry is known to operate for the different flavors in the quark- 
lepton system. In the context of the seesaw mechanism, such degeneracy is awkward to 
accommodate, at best. Furthermore, if and v T are nearly mass degenerate with a common 
mass of > 1/30 eV, then there is no simple way to address the solar neutrino puzzle through 
v e — or v e — v T oscillation. The possibilities appear bizarre - either v e is nearly degenerate 
with and z/ T , or there is a fourth "sterile" neutrino with the right mass and mixing 
parameters. In light of all this, it seems reasonable to consider alternatives sympathetically. 

We will argue now that large — v T oscillation angle can in fact arise, without requiring 
the neutrinos to be nearly degenerate, along the lines of some not entirely unsuccessful 
attempts to relate mixing angles and hierarchical masses in the quark sector. 

In the two sections following this one we shall analyze a complete model for the quark 
and lepton masses that is both overconstrained and phenomenologically successful (though, 
to be sure, its construction involves considerable guesswork). Since that analysis becomes 
rather complicated and perhaps intimidating, we shall first, in this section, highlight some 
salient features in a caricature model. 

A leading idea in many attempts to understand large mixing in the quark sector is 
to utilize the properties of special matrices, whose form might be constrained by simple 
symmetry requirements (e.g., requiring symmetry or antisymmetry) and selection rules. 
Such constraints can easily arise, as we shall see, from the group theory of unification, given 
specific choices for the Higgs fields whose condensation generates the masses. 

The first attempts along these lines, which remain quite intriguing and instructive, uti- 
lized a 2 x 2 system of the form 



Mi 



(° I) , l = M (5) 



with cif <C bf. Symmetric matrices arise, for example, in 5*0(10) (or left-right symmetric 
models |21[) where the matter fields are 16s if the relevant Higgs fields are 10s. The vanishing 
of the (1,1) entry can be ensured by a suitable flavor symmetry that distinguishes the two 
relevant families. The eigenvalues of this matrix are (mi, 777-2)/ — i a2 /b,b)f. Thus the off- 
diagonal element a/ is nearly the geometrical mean of the two eigenvalues. For convenience 
in discussion, we will refer to mass matrices of this form as type A. For type A mass matrices, 
the mixing angles in each sector, up and down, are given by 



tan# 



777! 



777 2 



(6) 

/ 



To obtain the observable mixing angle one must of course combine the mixing angles of the 
up and the down sectors. In a world with two flavors, this leads to the well-known formula 
for the Cabibbo angle P 



->c 





lm u 


V m s 


V m c 



(7) 



Using \Jrnd/rn s ~ 0.22 and ^m u /m c ~ 0.06, we see that Eq. (7) works within 30% for any 
value of the phase 0, and perfectly for a value of the phase parameter around ir/2. 
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A notable feature of the type A pattern, which remains valid even if the requirement of 
symmetry is relaxed, is that it can support a strong hierarchy of eigenvalues by means of a 
much weaker hierarchy of matrix elements. For example, with (a/b)f = 1/10, one obtains a 
large hierarchy {m\/rri2)f ~ 1/100. 

Now let us consider the implications of adopting the type A pattern for the [i — r sector, 
including the 2x2 Dirac mass matrices of the charged leptons (fi — r) and the neutrinos 
(u^ — Ur), and the Majorana mass matrix of the right-handed neutrinos (i/^ — i/^). Including 
the first family generally will not much affect the discussion of the \x — r sector, as we shall 
see. 

The three matrices in the \x — r leptonic sector have the form 

darY'" H/fR f° V 



with the understanding that in each case the (rr) entry dominates. One finds easily that 
the physical mass matrix for the light left-handed neutrinos, M^ 8ht = —M®(M^ I )~ 1 (M^ ) ) T , 
takes a similar form, viz. 

M iight _ 1 ( -dl T y V ( ) 

M » M R y*\-dl T y {dl T - 2d TT d, T y} J ' [J) 

Denoting the two eigenvalues as m V2 and m U3 , this yields: 



V 2 .._nu \ m v2 I™. 



s*v> tan ^ = \ ^T' ta < = ^- (io) 



Observe that the square-root formula holds for the mixing angles in all sectors, including 
the light neutrinos. 

To illustrate the possibilities, let us consider the phenomenologically interesting hierarchy 
ratio m Vti jm Vr ~ 1/10. With this ratio, one can accommodate m„ T ~ .03 eV as suggested by 
atmospheric neutrino oscillations and a value of m v consistent with the small-angle MSW 
solution (with very small m„J This ratio is achieved for (d^ T /d" T , y) = (1/5, 1/13). 
Combining the contributions from the charged lepton sector and from the neutrino sector, 
the physical oscillation angle following from Eq. (10) is 



nose 

V 



m-r 




(11) 



Although the mixing angle is not large in either sector (9^ T ~ ^Jm^/m T ~ 0.25 ~ 14° and 

9" lT ~ y/m^ /m„ T ~ yl/10 ~ 0.31 ~ 18°), if rj ~ 7T one obtains a near-maximal value for 
the physical mixing parameter. Actually the small angle approximation used in Eq. (11) is 
too crude; the precise expression is given by 

V V m r V m VT J V 



8111 29 » T = / , m A 2/ ^—T2 > (I 2 ) 

(1 + pi ) (1 + p^ 



where the ± corresponds to r) — (ir, 0). For 1] = it and m v ^jm VT = 1/10, sin 2 20°^ c ~ 0.79. 
This example shows that large oscillation angles can arise in a simple hierarchical model, 
without any extreme adjustment of parameters. As shown below, even larger oscillation 
angles can arise quite natually in more realistic models which can account also for V^,. 
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IV. A MORE AMBITIOUS 50(10) MODEL 







\v ch \ ~ 






V m t pi 



The model discussed in the previous section is not adequate for a unified treatment of 
the masses and mixings of the second and the third family fermions. Indeed, the square-root 
formula for V c b reads 

i /=r- /=n 

so that with ^Jm s /mb — 0.17 and ^Jm c /mt — 0.06, one cannot obtain observed value 
V^fe ~ 0.04 ± 0.003 for any value of x- Thus the simplest symmetrical type A mass matrices 
(Eq. (5)) cannot adequately describe the hierarchical masses and the mixings of the quarks. 

We now propose to study in detail a concrete proposal for asymmetric type A mass 
matrices that can be obtained within 5*0(10), and predicts correlations among the quark- 
lepton and up-down mass matrices that are phenomenologically acceptable. Before plunging 
into the analysis, let us briefly summarize its main conclusions. In this section, we shall 
temporarily ignore the first family. 

The four 2x2 Dirac mass matrices (U, D, L, and N) will be generated via four 
Yukawa couplings. The Majorana matrix of the right-handed neutrinos involves two 
additional Yukawa couplings. Thus there are seven parameters (six Yukawa couplings 
plus one ratio of vacuum expectation values) to describe 10 observables (eight masses - 
(m c , m t , m s , m b , m^, m T , m U2 ,m U3 ) - and two mixing angles - 9° s f and V c b) for the \i and the 
t families. The system is overconstrained, and predicts three relations among observables. 
Two of these concern the charged fermion sector, and they are reasonably well satisfied. 
The third prediction concerns the neutrino sector: the neutrino oscillation angle will be 
predicted as a function of the mass ratio m U2 /m u , A . For m U2 /m U3 = (1/10 — 1/30), consistent 
with mass determinations from the SuperK atmospheric neutrino data and the small angle 
MSW solution of the solar neutrino puzzle, a large [y^ — v T ) oscillation angle is obtained, as 
indicated by the SuperK atmospheric data. 

In forming hypotheses for the form of the couplings responsible for the mass matri- 
ces, we are guided by several considerations. We assume an underlying 50(10) unification 
symmetry, and that unified symmetry breaking occurs through a minimal system of low di- 



mensional Higgs multiplets - specifically, (45h), one pair of (16h) and ^16h/- A single 10h 
is employed for electroweak symmetry breaking. The hierarchical masses will be assumed 
to arise from type A mass matrices, as in Eq. (5), generalized to asymmetric form. In addi- 
tion, we must of course satisfy the broad phenomenological constraints m&(Mu) — m T (M\j), 
m s (Mu) 7^ ra M (Mu), and Vckm ^ 1- Combining these considerations with a restriction 
to low dimensional operators (d < 5), we are led to suggest the following set of Yukawa 



couplings [p3fl : 



^Yukawa = ^33 163 163 10 H + ^16 2 16 3 10 H 45 H + 

^16 2 16 3 16 h 16h + /123I62I63IOH • (14) 

Here M, a scale associated with the effective non-renormalizable interactions, could plau- 
sibly lie somewhere between the unification scale Mu and Mpi anc k. For example, the a 2 3 
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term might result from integrating out a 16 + 16 superfield pair with mass of order M; or 
M might be identified as Mpi anc k itself if the nonrenormalizable interactions are associated 
with gravity. A mass matrix of type A results if the first term, h 33 (10h), is dominant. This 
ensures m b (Mu) ~ m T (Mu) and m t (M v ) ~ mf (M v ). The remaining terms, responsible for 
off-diagonal mixings, must be smaller by about one order of magnitude. 

In more detail, our rationale for favoring the particular off-diagonal couplings 023 and #23 
is the following. If the fermions acquired mass only through 10s, one would obtain too much 
symmetry between down quark and charged lepton masses, and in particular the bad rela- 
tion m s (Mu) = m /J (Mu). To avoid this, while eschewing proliferation of extraneous Higgs 
multiplets, the simplest possibility is to bring in the vacuum expectation value of the (1,1,15) 
component (under G224) in the 45h that is proportional to {B — L). A vacuum expectation 
value of this form figures prominently in the Dimpoulos-Wilczek mechanism for doublet- 
triplet splitting in SO (10) H, as we shall discuss at length in connection with proton decay. 



Now if one restricts to d = 4 and d = 5 operators in the superpotential, the only relevant 
5*0(10) invariant effective coupling that involves the 45h is our ajjl6il6jlOH45H,M7. In 
this term, only the 120 in the decomposition 10h x 45h (= 10 + 120 + 320) can contribute 
to fermion mass matrix, and its contribution is antisymmetric in These two couplings 

still do not distinguish up and down quark mass patterns, and so if they were the whole 
story one would have a trivial, identity, Cabibbo-Kobayashi-Maskawa (CKM) matrix. The 
effective couplings gjjl6il6jl6Hl6H/7V7 remedy this problem, without requiring addition to 
our small set of Higgs multiplets. When a vacuum expectation value of order the unification 
scale for the £77(3) x £77(2) x (7(1) singlet component of 16h and a vacuum expectation 
value of order the electroweak scale for the SU(2) x £7(1) breaking component of 16h are 
inserted, this term contributes to quark and lepton mass matrices. To be more precise, since 
the electroweak doublet contained in 16h has the quantum numbers of a down-type Higgs 
doublet, it contributes to the down quark and charged lepton mass matrices, but not to the 
up sector. The resulting up-down asymmetry generates non-zero CKM mixing angles. 

With these four effective Yukawa couplings, the Dirac mass matrices of quarks and 
leptons of the second and the third families at the unification scale take the form: 

tt ( e + a \ n ( e + r A 

U= {-e + a 1 ) mu > D= [-e + V 1 ) m ^ 

at ( "3e + <x\ T ( -3e + rA 

N= W + a 1 ) m ^ L= {3e + V 1 ) mD - (15) 

Here the matrices are multiplied by left-handed fermion fields from the left and by anti- 
fermion fields from the right. (U, D) stand for the mass matrices of up and down quarks, 
while (TV, L) are the Dirac mass matrices of the neutrinos and the charged leptons. 

The entries (l,e, a) arise respectively from the ^33,023 and terms in Eq. (14), while 
r) entering into D and L receives contributions from both (723 and ^23; thus 77 ^ a. Note the 
quark-lepton correlations between (U,N) as well as (D,L), and the up-down correlation 
between (U, D) as well as (TV, L). These correlations arise because of the symmetry structure 
of £O(10). The relative factor of —3 between quarks and leptons involving the e entry reflects 
the fact that (45h) oc (B — L), while the antisymmetry in this entry arises from the £O(10) 
structure as explained above. 

Assuming e,rj,a 1, we obtain at the unification scale: 
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m c 


1 2 2| 

— e — cr | , 


m t 




rrifj, 


~ |9e 2 — 77 2 1 


m T 


\v cb 


~ |cr — T]\ . 



iH^i ~ i e 2 _ „2i 



lm-.Hl - 8e" 



(16) 



One virtue of the asymmetric nature of the mass matrices (Eq. (15)) is worth noting. The 

simple expression | cr — 77 | for | V c b | actually stands for ^m s /rrib (jf~) — ^Jm c /m t 

Comparing with Eq. (13), which was based on symmetric type A mass matrices and which 
led to too big a value of \V c b\, we see that the square-root mass ratios are multiplied by the 
respective asymmetry factors (^zf) 1 ^ 2 and (frf) 1 ^ 2 (see Eq. (15)). For e relatively negative 
compared to rj and a (as suggested by considerations of the leptonic mixings, see below), 
these factors provide the desired lowering of V c b compared to Eq. (13). 

In Eq. (16) all the mass and mixing angle parameters are to be identified with those at 
the unification scale. One can evaluate a, r], e in terms of the observed masses and mixing 
(extrapolated to the unification scale): 



° — l v cb + m s /m b - m c /m t ]/(2V cb ) , 
V — \~ V cb + m s/m b - m c /m t }/ (2V cb ) , 
; 2 ~ {a 2 + m c /m t } ~ {rj 1 + m A1 /m r )/9 



This leads to the sum rule: 



m t 



5- 
4 



cb 



+ l^v 

16 cb 2 m T 



9 m L 
2 m t 



.1/2 



(17) 



(18) 



A word of explanation is needed here. The parameters a, 7] are in general complex. In 
Eqs. (16)-(18) therefore, one should interpret the mass ratios and V c b to be also complex. Let 
us denote m c /m t = r] ct \m c /m t \,m s /mb = i] s b\m s /m b \,V c b = r] cb \V cb \, with |^| = 1. Noting 
that each term on the RHS of Eq. (18) is small compared to m s /mb — (1/30 — 1/50) (e.g: 
m c /m t ~ 1/300, iVcbWm^/ (2m T )\ 1 ^ 2 ~ 1/140), we see that essentially only one choice of the 
phase factors can possibly allow the sum rule to work so that the RHS is maximized - i.e., 
r/ ct = —l,r)n T = +1 and r] sb = —1 with {±V^,[ J 1 / 2 } having a negative sign. Thus, all phases 
are constrained to be near the CP conserving limit (0 or it). For simplicity, we will take the 
rjij's to be real. As long as the phases are small (< 10°, say), they would not of course affect 
our results. 

The relatively simple pattern shown in Eq. (15) provides a reasonable fit to all the masses 
and the mixing of the quarks and the charged leptons in the second and the third families. 
For example, if we take as input mf hyslcal = 174 GeV, m c (m c ) = 1.37 GeV, V c b = 0.045, in 
agreement with the values advocated in Ref. 



25|j26[] , and the known and r lepton masses, 



then we obtain the predictions 



mb(mb) ^ 4.9 GeV , 
m s (l GeV) ~ 116 MeV . 



(19) 



In quoting the numbers in Eq. (19), we have extrapolated the GUT scale values down to low 
energies using the beta functions of the minimal supersymmetric extension of the Standard 
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Model (MSSM), assuming a s (M z ) = 0.118, an effective SUSY threshold of 500 GeV and 
tan/3 = 5. Our results depend only weakly on these input choices, so long as tan/3 is neither 
too large (> 30) nor too small (< 2). 

The mass of the strange quark is somewhat low compared to the central value advocated 
in Ref. ||25|| , but is closer in value to the recent lattice determinations P7| . The fe-quark 



mass prediction is also in reasonable agreement with determination from T spectroscopy. In 
this regard, it should be mentioned that exact b — r unification in supersymmetric unified 
models, without taking finite threshold effects due to the gluino into account, would yield 
a 6-quark mass which is about 10-20% above the experimental value for a wide range of 
the parameter tan/3 |[28|| . In our case, b and r masses are not exactly equal at Mu, m b is 
about 8% lower than m T (see Eq. (16)). This difference, which has its origin in the (B — L) 
generator associated with the off-diagonal entry e in Eq. (15), leads to better agreement 
with the experimental value of m^m;,). 

The parameters a, rj, e, mu and m D are found to be 

a ~ —0.110 r] cb , 7] ~ —0.151 rj cb , e ~ 0.095 m , 
mu ~ m t {M v ) ~ (100 - 120) GeV , m D ~ m b {M v ) ~ 1.5 GeV . (20) 

(m = ±1 is the phase of e.) In addition to the two predictions in Eq. (19), the Dirac masses 
of the neutrinos and the left-handed mixing angles in the fi — r sector are determined to be 

m£(Mu) ~ m t (M v ) ~ 100 - 120 GeV , 
m^(Mu) ~ (9e 2 - a 2 )mu ~ 8 GeV , 

6l T ~ -3e + ij ~ -0.437r/ £ (for Vch /r] e = +1) . (21) 

Vcb/Ve = +1 i s required to maximize 9 . Note that the Dirac mass of is quite different 
from m c (My) ~ 300 MeV. This is contrary to what is often adopted in 5O(10). The 
difference arises because of the type A pattern of the mass matrices and the factor of —3 
associated with the (B — L) generator that goes into the Dirac mass [25 . 



Given the bizarre pattern of quark and lepton masses and mixing, we regard the overall 
fit to all of them, good to within 10%, using the pattern shown in Eq. (15), as reason to 
take this pattern seriously. 

Note that owing to the asymmetric mass matrix, the square-root formula for the mixing 
angle 9 £ receives a correction given by the factor [(— 3e + r))/(3e + r])} 1 ^ 2 - For the values of 
rj and e determined in Eq. (21), this factor is about 1.8 for rj e /rj c b = +1. 

Although all the entries for the Dirac mass matrix are now fixed, to obtain the parameters 
for the light neutrinos one needs to specify the Majorana mass matrix. Unfortunately, here 
there is much less information to guide our hypotheses. For concreteness, let us imagine 
that this too takes the type A form: 

M ?=( y \) M ^ ( 22 ) 

where we allow y = i] y \y\ to have either sign, i.e., rj y = ±1. Note that Majorana mass 
matrices are constrained to be symmetric by Lorentz invariance. The seesaw mass matrix 
(—N(M^)~ 1 N T ) for the light (i/„ - v T ) system is then 
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MT=(° A i)^, (23) 



A BJ M R ' 

where A ~ (a 2 - 9e 2 )/y and B ~ -(a + 3e)(cx + 3e - 2y)/y 2 . With A < B, this yields 



M/j m u , B 2 ^ V 



Correspondingly, 



A/ rrt„„ \ I I 



3|e| + lal ±2,/^ 



(25) 



where ?7 e = d=l and we have used the fact that e and a are relatively negative for 
Vcb/Ve — +1 (See Eq. (21)). For a given choice of the sign of y relative to that 
of e, and for a given mass ratio m U2 /m V3 , we can now determine yr] e using Eq. (25) 
and the values of e and a obtained in Eq. (20). Corresponding to m U2 /m U3 = 
(1/10,1/15,1/20,1/25,1/30), we calculate (yr] e )+ = (0.0543,0.0500,0.0468,0.0444,0.0424) 
and (yr] e )- = (0.2359,0.3781,0.7681, 8.4095,-1.0519), where the subscripts ± correspond 
to r]yr] e — ±1. The case of rj y rj t = +1 typically requires smaller values of \y\ than for the case 
of r]yr] e — — 1. The former (r) y r) e = +1) is more in accord with the idea of the type A matrix 
with flavor symmetries being the origin of hierarchical masses, than the latter (r) y r) e = — 1). 
We obtain for the neutrino oscillation angle: 



V 0.437 ± J-^. (26) 

If r)y/r) e = +1, the relative + sign should be chosen, and for r) y /r) e = —1 the relative — sign 
should be chosen. Taking the relative + sign, i.e., r) y r) t = +1 (in accord with the smaller 
values of y < 1/10), and using the more precise expression given in Eq. (12) (and replacing 
\Jm^/m T by the numerical value of the charged lepton mixing angle ~ 0.437), we obtain 

sin 2 2^ T = (0.96,0.91,0.86,0.83,0.81) 
for m v Jm v . A = (1/10, 1/15, 1/20, 1/25, 1/30) . (27) 

As previously advertised, we see that one can derive rather plausibly a large — v T 
oscillation angle sin 2 2# MT > 0.8, together with an understanding of hierarchical masses and 
mixing of the quarks and the charged leptons, while maintaining a large hierarchy in the 
seesaw derived masses (m U2 /m U3 = 1/10 — 1/30) of and v T , all within a unified framework 
including both quarks and leptons. In the example exhibited here, the mixing angles for 
the mass eigenstates of neither the neutrinos nor the charged leptons are really large, ~ 
0.437 ~ 23° and 0£ T ~ (0.18 - 0.31) w (10 - 18)°, yet the oscillation angle obtained by 
combining the two is near-maximal. This contrasts with most previous work, in which 
a large oscillation angle is obtained either entirely from the neutrino sector (with nearly 
degenerate neutrinos) or entirely from the charged lepton sector. In our case, the mass 
eigenstates of the neutrinos and the charged leptons are approximately also the respective 
gauge eigenstates. 
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V. INCLUSION OF THE FIRST FAMILY: u E - OSCILLATION 

There are several alternative ways to include the first family, without upsetting the 
successful predictions of the 2-3 sector. In the absence of a deeper understanding, the 
theoretical uncertainties in analyzing the masses and mixings of the first family are much 
greater than for the heavier families, simply because the masses of the first family are so 
small, that relatively tiny perturbations can significantly affect their values. With this 
warning, we will now briefly consider, as an illustrative example and "proof of principle" , a 
minimal extension to the first family, inspired by the type A pattern, and also a noteworthy 
variant. 

The 3x3 Dirac mass matrices in the four sectors take the form: 

e' + rj' \ 

e + r] \ m D , 
-e + r] 1 / 

Se' + r]' \ 
3e' + rf -3e + r? m D . (28) 

3e + r] I J 

At the level of underlying primary couplings, the e' term is very similar to the e term, 
arising through a Yukawa coupling ai2l6il6245Hl0H- The rj term arises through the cou- 
pling (7i 2 16i16216h16h- With e,a,r],mu and mo determined essentially by considerations 
of the second and the third families (Eq. (20)), we now have just two new parameters in 
Eq. (28) - i.e. e' and rj - which describe five new observables in the quark and charged 
lepton sector: m u , m d , m e , 6c and V ub . Thus with m u « 1.5 MeV (at Mu) and m e /m^ taken 
as inputs to fix e' and r/, one can calculate the other three observables. In addition, 
will also be determined. 

At the outset, it is worth noting that this specific pattern (Eq. (28)) is consistent with 
the empirical unification scale relations: ~ 3m e ,m s ~ [[50] and m& ~ m T , which 

in turn imply m d m s m b ps m e m^m T1 and thus Det(D) ~ Det(L). Eq. (28) obeys this 
relation to a good approximation, for the following reason. From Det(U) = e' 2 mfj, one 

obtains e' ~ y 1 m u /m c {m c /m t ) ~ 2 x 10 -4 . Using Det(L) = (9e /2 — r]' 2 )m 3 D = mem^mr, 

and neglecting the e' 2 term (justified a posteriori), one obtains \rf\ ~ \Jm e /m^{m^/m T ) 
~ 4.4 x 10 -3 . Since |e'| <ti \rj'\, one gets Det(D) ~ Det(L), as described. 

Combining the two predictions for the second and the third families given by Eq. (19), 
we are thus led to a total of five predictions for observable parameters of the quark and 
charged lepton system of the three families. 

m^mft) ~ 4.9 GeV 
m s (l GeV) ~ 116 MeV 
m d (l GeV) ~ 8 MeV 

6c- 



\V ub /V cb \ ~ ^m u /m c ~ 0.07 . (29) 
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Further, the Dirac masses and mixings of the neutrinos and the mixings of the charged 
leptons also get determined. Including those for the /i — r families listed in Eq. (21) we 
obtain: 

m£(Mu) ~ 100 - 120 GeV; mJ(M D ) ~ 8 GeV, 6^ ~ -.0.437^ 
m° ~ [9e /2 /(9e 2 - a 2 )\m v ~ 0.4 MeV 



0' ~ 



77 — 3e 



\Jm e /m^ ~ 0.85^/me/m^ ~ 0.06 



C ~ ^meKK/ra,) ~ 0.0012 . (30) 

In evaluating we have assumed e' and 77' to be relatively positive. 

Note that the first five predictions in Eq. (29) pertaining to observed parameters in the 
quark system are fairly successful. Considering the bizarre pattern of the masses and mixings 
of the fermions in the three families (recall comments on V c b, my ) jm T ^ m s /m^ and m^/m e ), 
we feel that the success of the mass pattern exhibited by Eq. (28) is rather remarkable. It 
is worth stressing that (a) the asymmetric mass pattern induced by the 023 coupling (see 
Eq. (14)), (b) its contribution being proportional to B — L and (c) the g 2 3 coupling leading 
to non-trivial CKM matrix have all played essential roles in leading to a successful mass 
pattern with a minimal Higgs system - i.e., involving a single 45#, 16h and 10h - two of 
which are needed for unification scale symmetry breaking of 5*0(10) anyhow. This is one 
reason for taking patterns like Eq. (28) seriously as a guide for considerations on proton 
decay. A particularly interesting variant is obtained in the limit e' ^ 0, as we will discuss 
at the end of this section. 

Turning to the Majorana mass matrix, one could in general introduce (1,1), (1,2) and 
(1,3) entries involving the first family. We choose the (1,2) element to be zero, both because 
by itself (i.e., with (1,2) 7^ 0, but (1,1) = (1,3) = 0) it would upset the success of the 

— v T sector obtained above, and also for economy in parameters. We therefore consider 
the following pattern: 



M, 



R 



fx z\ 

y 
\z y 1/ 



M R . (31) 



Consistent with our presumption that hierarchical masses have their origin in the type A 
pattern, we will assume that x ~ z 2 and z < y < 1 and that the zero in the (2,2) entry only 
reflects that this entry is small compared to y 2 . 

Note that all entries in the Majorana matrix arise through /jjl6jl6jl6Hl6H, which 
however do not contribute to the Dirac mass matrices. As we have stressed before, they 
therefore need not have the same flavor structure as the Dirac mass matrices of Eq. (28). 
Given the Dirac matrix N, we can work out the light neutrino mass matrix after the seesaw 
diagonalization. It is given by: 

M^ ht = (32) 
M R xy 2 v ; 
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/ 9e' 2 (x-z 2 ) 3e'y{-3e'z + (-3e + a)x} 3e'{xy - (3e + a)(x - z 2 )} 

3e'y{-3e'z + (-3e + a)x} -9e' 2 y 2 (3e + a)y{3e'z - (-3e + a)x} 

\3e'{xy - (3e + a)(x - z 2 )} (3e + a)y{3e'z - (-3e + (3e + a){-2xy + (3e + - 2 2 )} 

Three features arising from Eq. (30) are especially noteworthy. First, that the parameter 
x figures prominently only in the mass of z/ e , which is given by 

81e' 4 ,771* KT 10 eV /QQ , 
-) ~ , (33) 



e xy 4 (9e 2 -a 2 ) 2X M R J x 

where we have inserted the mass of v T ~ 1/20 eV. For x = (10~ 5 — 10~ 4 ), m Ue ps (2 x 10~ 5 — 
2 x 10~ 6 ) eV. (x should be less than 10~ 4 , otherwise one would run into a conflict with 
proton lifetime, see later.) Second, the discussion of the \i — t sector in Sec. IV is corrected 
by terms that are proportional to e', which is tiny, or by terms proportional to z. As long 
as z 2 < x/3 (from the (3,3) entry) and z < 10 3 x (from the (2,3) entry), the discussion 
of the (2, 3) sector will remain essentially unaltered. And third, that the v e — mixing 
angle arising from the light neutrino mass matrix itself (as opposed to the contribution from 
charged lepton mixing) is tiny, but still significant for the small angle MSW explanation of 
the solar neutrino puzzle. The v e — v T mixing is also small, since e' ~ 2 x 10~ 4 is extremely 
small. 

Note that the 3x3 Majorana mass matrix Mjf introduces four new parameters: x,y,z 
and M R , of which the magnitude of M R ~ (few to 10) x 10 14 GeV can quite plausibly be 
justified in the context of supersymmetric unification (see Sec. II). There are altogether six 
observables in the three neutrino system: the three masses and the three oscillation angles. 
While these depend on parameters of both the Majorana and the the Dirac mass matrices, 
since the latter are fully determined from the quark and the charged lepton system, one can 
expect at least two predictions for the neutrino system. These can be taken to be 9°^ and, 
for example, m Ue . In addition, to the extent that the magnitude of M R can plausibly be 
justified, that of m Vr ~ (1/10 — 1/30) eV can also be anticipated. We should add that for 
the limiting case of e' = (a parameter of the Dirac sector), one would obtain predictions 
for 9°^ and 9™° as well, regardless of the values of x, y, z (see below). 

As an example, let us take x = 10~ 4 ,y = 0.05 and z = —0.002 along with the "central 
values" of a, r],e given in Eq. (21) and e' = 2 x 10 -4 ,?/ = 4.4 x 10~ 3 (from a fit to the 
■u-quark and electron masses). Let us also put, in accord with discussions in Sec. II, 
M R Ki (5 — 15) x 10 14 GeV. The mass eigenvalues and the mixing angles in the neutrino 
sector are then: 



m Ur ^ (1/10 - 1/30) eV 

~ 6.7 x 10~ 3 (1 - 1/2) eV 
m Ue ~2x 10~ 6 (1 - 1/2) eV 

9^ T ~ 0.25, 9 v eil ~ 0.009, 9 v eT ~ 0.0034 . (34) 

Combining with the mixing angles in the charged lepton sector, the oscillation angles, in- 
cluding that for z/ M — u T sector given before, are given by (setting the relative phases to be 
or 7r): 
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6^ ~ 0.437 + ^m v Jm v 

6{ T - 0£ T ~ 5.5 x 10" 4 ± 0.0034 . (35) 



"v 3 



Though the mixing angles involving the first family are small, Q v does play an important 
role for solar neutrinos. Combining Q v with the corresponding mixing angle from the charged 

lepton sector, 9*^ ~ 0. 85 ^rrie/m^ ~ 0.06, where the factor 0.85 is a correction factor from 
the asymmetric e' term in Eq. (28), we find that 9° 8 ° ~ 0.05 if the two contributions subtract. 
This is just about the right value to explain the solar neutrino data via the small angle MSW 
mechanism Note that without the (1,3) entry z in Eq. (31), Q v would have been less 
than 0.002, which would have resulted in 6° s ^ being larger than the required value for MSW 
by about 15-20%. This is a possible motivation to introduce the (1,3) entry z into Eq. (31). 
Since the discussion of the (2-3) sector is little affected so long as z < 10 3 x and z 2 < x/3, 
we will allow a range \z\ = (0.002 — 0.006), corresponding to x — (10~ 5 — 10~ 4 ). This will 
be significant for the discussion of proton decay, especially for final states involving charged 
leptons. 

Another consideration for the choice of the various entries in is the generation of 
cosmological baryon asymmetry. Note that for x = 10~ 5 — 10 -4 , the mass of is in the 
range 5 x (10 9 — 10 11 ) GeV, which is of the right magnitude for producing following 
reheating and inducing lepton asymmetry in decay, that is subsequently converted to 
baryon asymmetry via the electroweak sphalerons For the lepton number violating 

decay of z/f to be out of equilibrium, it is necessary that the mixing of with be less 
than a few times 10~ 3 . Otherwise the decay — > H°+u^, which can proceed by utilizing the 
large (order one) Dirac Yukawa coupling of will bring uf into equilibrium, inhibiting 
the generation of lepton asymmetry. The choice z ~ 0.002 is consistent with this constraint. 

While we cannot claim that our several choices have been unique, we have demonstrated 
that a rather simple pattern for the four Dirac mass matrices, motivated and constrained 
by the group structure of 5O(10), is consistent within 10 to 20% with all observed masses 
and mixing of the quarks and the charged leptons. This fit is significantly overconstrained, 
as we have discussed. The same pattern, supplemented with a similar structure for the 
Majorana mass matrix, quite plausibly accommodates both the SuperKamiokande result 
with the large — v T oscillation angle required for the atmospheric neutrinos and a small 
— V/j. oscillation angle relevant for theories of the solar neutrino deficit. 

The e' = Limit: A special case 

Before turning to proton decay, it is worth noting that much of our discussion of fermion 
masses and mixings, including those of the neutrinos, is essentially unaltered if we go to the 
limit e' — > of Eq. (28). This limit clearly involves: 



m u = 0, 6 C ^ \Jm d /m s 



11 — e -Jm d /m b {m s /m h ) ~ (2.1)(0.039)(0.023) ~ 0.0019 
V V + e v 

0, &L = ei T = . (36) 
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All other predictions will remain unaltered. Now, among the observed quantities in the list 
above, 9c — \Jmd/m s is indeed a good result. Considering that m u /m t ~ 10" 5 , m u = 
is also a pretty good result. There are of course, plausible small corrections arising from 
higher dimensional operators (for example) involving Planck scale physics which could induce 
a small value for m u . The question arises: Does such a small correction contributing to m u 
arise primarily through the (1,2) entry (i.e., e' 7^ 0) in U and N or does it arise primarily 
through the (1,1) entry - to be called 5 7^ if e' = - of U and A? In case of the former (with 
5 = 0), one would get e'~2x 10~ 4 , as noted above; while for the latter (with e' = 0, 5 7^ 0), 
one would obtain 5 ~ m u /m t ~ 10 -5 . Note that putting e' = does not matter for D and L 
because rf ^> e'. The two cases therefore yield essentially identical results as regard fermion 
masses and mixings. Differences of about 20% between their predictions can arise only as 
regards 9° B ^ for which the former can fare better (depending upon the choice of phases) 
than the latter; while the latter (5 7^ 0, e' = 0) fares better than the former as regards 
9c- One possible way to distinguish these two closely related variants for fermion masses 
and mixings will be a more precise determination of V ub . The former predicts it to be near 
^m u /m c V cb ~ (0.07) (0.04) ~ 0.0028, while the latter is somewhat smaller (V ub ~ 0.0019). 
We will refer to these two variants as cases I and II respectively. 

Case I : e' « 2 x 10" 4 , 5 = 

Case II : 5 « 10~ 5 , e' = (37) 

As we will see, however, proton decay amplitudes differ significantly between the two 
cases. Typically they are larger for case I (by a factor of 2-3) than for case II. Not knowing 
the precise origin of the small corrections contributing to m u , we will consider cases I and 
II as equally viable guides for forming expectations for proton decay. We will in fact derive 
the proton decay amplitudes with e' 7^ (case I). Given that 5 <C e', we obtain the results 
for case II simply by setting e' ~ 0. 



VI. NEUTRINO MASSES AND PROTON DECAY 



In an earlier paper [12] we pointed out that the theory of neutrino masses can significantly 
affect expectations for proton decay, as regards both its rate and, especially, its branching 
ratios. This happens because in supersymmetric unified theories a new set of color triplet 
fields is needed to generate heavy Majorana masses for the right-handed (RH) neutrinos, as 
required for the seesaw mechanism. Exchange of its superpartners generates new dimension 
5 operators, that appear in addition to the "standard" d = 5 operators in the effective 
Lagrangian. The standard operators arise through the exchange of color triplets that are 
related to the electroweak doublets, such as those appearing in the 5 + 5 of SU (5) or the 

10 of so (10) hhhh. 



The standard d = 5 operators are estimated to dominate over the gauge boson mediated 
d = 6 proton decay operators. Furthermore, owing to a combination of color antisymmetry, 
Bose symmetry of the superfields, and hierarchical Yukawa couplings of the fermions, they 
predict that the dominant modes are V ll K + and T^vr 4 " modes, while e + 7r° and e + K° and 
even /i + 7r° and fi + K° modes are highly suppressed, at least for small and moderate values 



of tan/? (< 15) JHM-P 
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In Ref. |T|] we argued that for neutrino masses in a plausible range the new d = 5 op- 



erators could compete favorably with the standard ones, and induce potentially observable, 
though not yet excluded, proton decay rates. Within the new contribution charged lepton 
decay mode amplitudes, including those for fi + K°, /i + 7r° and possibly e + K° and e + n°, were 
generally found to be significant relative to the VK + and Vtt + modes, even for low tan/5. 
Thus charged lepton modes of nucleon decay potentially furnish an independent signature 
for the interactions that play a central role in the theory of neutrino masses. 

As mentioned earlier, in our previous quantitative work we had to hypothesize values for 
the neutrino masses. Now, after the SuperK observations, we can revisit the quantitative 
aspect with some crucial experimental information reliably in hand. Changes are induced 
not only for the new d — 5 operators that are directly related to neutrino masses, but also 
for the standard d = 5 operators. 

To address these issues in detail, we must specify the nature of the color triplet Higgsino 
couplings. This issue is closely tied to the mechanism of doublet -triplet splitting in 5*0(10), 
to which we now turn. 



A. Natural doublet-triplet splitting in 50(10) 

In supersymmetric 50(10), a natural doublet-triplet splitting can be achieved by cou- 
pling the adjoint Higgs 45h to a 10h and a 10' H , with 45h acquiring a unification-scale 



VEV in the B — L direction f24]| : (45h) = (a, a, a, 0, 0) x r 2 with a ~ M\j. As discussed in 
Section II, to generate CKM mixing for fermions we require an (16h)^ that acquires an elec- 
troweak scale vacuum expectation value. To insure accurate gauge coupling unification, the 
effective low energy theory should not contain split multiplets beyond those of MSSM. Thus 
the MSSM Higgs doublets must be linear combinations of the SU{2) L doublets in 10h and 
16h- A simple set of superpotential terms that ensures this and incorporates doublet-triplet 
splitting is: 



W H = A10h45 h 10 , h + M 10 10' H 2 + A'16 h 16 h 10h + M 16 16 H 16 H • (38) 



A complete superpotential for 45h, 16h, 16h, 10//, 10' H and possibly other fields, which 
ensure that 45h, 16h and 16h acquire unification scale VEVs with (45h) being along the 
(B — L) direction, that exactly two Higgs doublets (H u , Ha) remain light, with Hd being a 
linear combination of (10n)d and (16n)d, and that there are no unwanted pseudoGoldstone 
bosons, can be constructed [p5H 39|| . [] The various possibilities generate different predictions 
for threshold corrections in the unification of gauge couplings, for example. As we will 
explain, such differences will have implications for proton decay rate, and we will allow for 
such effects. 

The Higgs doublet and the color triplet mass matrices following from Eq. (38) are, in 
SU(5) notation, 



4 It is intriguing that the Higgs fields used in [37 can be neatly incorporated into a single adjoint 
representation of E7. 



18 



no H 
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(39) 



With the vacuum expectation value (45h) in the B — L direction it does not contribute 
to the doublet matrix, so one pair of Higgs doublet remains light, while all triplets acquire 
unification scale masses. The light MSSM Higgs doublets are 



H u = 10 u , H d = cos7lO d + sin7l6 d 



(40) 



with tan7 = A' (16h) /Miq. Consequently, (10) d = 00S7 v d , (16d) = sury v d , with (Hd) = Vd 
and (16 d) and (10d) denoting the electroweak VEVs of those multiplets. Note that the H u 
is purely in IOh and that (lO^) 2 + (16 d ) 2 = v\. 

This pattern of gauge symmetry breaking, while motivated on separate grounds, nicely 
harmonizes with the fermion mass pattern advocated in Sees. IV and V. Specifically, the 
gij terms will contribute to the down-flavored fermion masses, while there are no analogous 
terms in the up-flavor sector. We also note the relation tan7 tan/3 ~ mt/m& ~ 60. Since 
tan (3 is separately observable, the angle 7, which will prove relevant for proton decay, is 
thereby anchored. 



B. Baryon number violation 

By combining Eqs. (38)-(40) with the Yukawa couplings from Eq. (14), we can now 
obtain the effective baryon number violating superpotential. Denote the color triplet and 

anti-triplet in IOh as (H C ,H C ), in 10' H as (H' C ,H' C ) and in (16h, 16h) as (H C ,H C )- 
The relevant Yukawa couplings of these color triplets, extracted from Eq. (14) aref] 
(i) KAlQiQsHc + Q l L J Hc), (h) -a^Q.QjH^B - L) Qj + Q % L^(B - L) L .], (iii) 

(g^ (16h) /M)QiLjH c , and (iv) ^16h ) /M )QiQjHc- Note that for the coupling 
there are two contractions: IOh x 45h D 10 + 120. Whereas the antisymmetric 120 con- 
tributes to fermion masses (for (45h) oc B — L), as in Eq. (15), it is the symmetric 10 
that leads to the Yukawa couplings of H c and He- Although the term can arise through 
quantum gravity or stringy effects, for concreteness, we have assumed that it arises by in- 
tegrating out a pair of superheavy 16 + 16 states which couple through the renormalizable 
interactions W D (16 2 16)10 H + (16 3 T6)45 H + Myl6T6.Q This is why the (B-L) generator 
appears in the couplings involving He, He- The strength of this coupling is then fixed 
in terms of the corresponding doublet coupling. Similarly, in the coupling, there are two 
5*0(10) contractions, and we shall assume both to have comparable strength. 

Baryon number violating operators of the type QQQL/M are the ones that connect to 
ordinary quarks and leptons by wino dressing, and which generally dominate (but see below). 



5 Here we have absorbed the factor (45#) /M into a^- and the factor (16 h) /M into gij. We use the 
same notation (ajj and g^) for these redefined quantities. As for fy, we define fij = fa (16h) /M. 

6 A variant is to interchange the indices 2 and 3 in this renormalizable interaction. 
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Integrating out the color triplet fields, one arrives at the following effective superpotential 
terms involving these operators: 

= M~q[(u t Hd') {u T HV'£ - d' T HV'u'} + 3(u T Hd') {u T AV'l - d! T AV'v') 

- (u T Ad') {u T HV'£ - d' T HV'u) - 3(u T Ad>) {u T AV£ - d' T AV'u'} 

- tan7( M T iM , ) {u T GV'£ - d ,T GV'u'} + tan 7 (« T 'Ad') {u T GV'£ - d' T GV'u'}} 

+ M^{u T Fd') {u T GV'l - d' T GV'u'} . (41) 

Here M cS = (Aa) 2 /M 10 . u and t denote the column matrices of the physical left-handed 
up quark and charged lepton superfields in the supersymmetric basis (i.e., a basis in which 
neutral gaugino interactions are flavor diagonal). The d' and v' fields are related to the 
physical down quark and light neutrino fields by the CKM matrices for quarks and leptons: 
d' = VcKud and v' = Vq KM u, while V = V^V^, where V u and diagonalize respectively the 
left-handed up quark and charged lepton mass matrices: = V u vS m \ where (g) and (m) 
denote the gauge and mass eigenstates. In writing Eq. (41), the color indices (a,/3, 7) on 
quark fields are suppressed and use is made of the fact that (ti^Hd'p — d'J Hup) = 
which holds because of antisymmetry under the interchange a <-> (3 and because H is 
symmetric. The superscript (L) on signifies that all the fields in belong to 

SU(2)l doublets. We comment later on the contributions from involving RRRR 

operators of the form u c u c d c e c , involving SU(2)l singlets, which can be important in certain 
range of supersymmetric parameter space. 

The 3x3 matrices (H, A, G, F) operate in the family space and are related to the Yukawa 
coupling matrices given by the elements (hij , , gij , fy) respectively, as follows: 



{H,A,G,F) = Vj(h,a,g,f)V v 



(42) 



The matrices (h, a, g, f) are related to those appearing in Eqs. (28), (31) and are given by: 



/o 0\ 



h = h. 




(43) 



where = j^v^jM. Values of a,r),e,r)',e' and y have been obtained in Sec. IV-V from 
considerations of fermion masses and mixings. From proton decay constraints and the 
neutrino sector we have estimated x ~ (10~ 5 — 10~ 4 ) and z ~ (0.002 — 0.006). V u can be 
worked out using Eq. (28) and (20), and is found to be close to an identity matrix. Its 
largest off-diagonal entry is the (1,2) element ~ e'/( e2 — °" 2 ) — —0.06. However the matrix 
Vi that diagonalizes the charged lepton mass matrix L is far from trivial, and so V = V^Ve 
picks up a substantial (2,3) element ~ — 3e + rj ~ —0.437, reflecting sizable fi — r mixing. 
The numerical values of the matrices V u , Vt and V are given in Appendix A, as are the 
matrices (//. .4. G. / '). 
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Even if one is skeptical of our particular pattern of fermion mass matrices, it is difficult 
to avoid the general conclusion^ that if neutrino masses are hierarchical, a large — v T 
oscillation angle (the SuperK result) requires a sizable fi — r mixing angle (say > 0.3). This 
has significant implications for proton decay, as we now discuss. 

After wino dressing, which dominates over gluino dressing for tan/5 < 20, each of the 
seven terms in Eq. (41) leads to twelve four fermion operators for proton decay into 77+ X, a 
subset of which was exhibited in Ref. |]12|| . In Appendix B.l we give the complete expression 
for the neutrino as well as the charged lepton decay modes of the proton. Representative 
contributions to the proton decay amplitudes are analyzed in more detail for the dominant 
modes in Appendix B.2. We now briefly summarize the results of a lengthy investigation of 
the net effect of all these terms. 

To evaluate the strength of each term in Eq. (41), we need h 33 and ^33 (see Eq. (43)) 
at M\j. h 33 is determined using h 33 v u ~ m t (Mu) — 100 — 120 GeV, which yields h 33 ~ 1/2. 
/33 can be determined as follows. Using Eq. (23) for the mass matrix of the light z/^ — v T 
sector, we have 

m V3 ~ Bml/Ma, (44) 

where 

B = -( a + 3e)(a + 3e-2y) ~ 5. (45) 

Here we have put a and e from Eq. (20) and used y = 0.047, corresponding to m V2 /m Uz = 
1/15. Putting 777,1,3 = (1/20 eV)C, where £ = 2 to 1/2, corresponding to SuperK results, 
and mu m t (M v ), we get M R ~ 10 15 GeV/C Using Mr = f 33 v 2 R /M (see Eq. (4)), with 
vr = (16 H ) = (2 x 10 16 ) GeVKR and k r « 1/2 to 2, we find 

/as = fssVR/M ~ (1/20)(1/«hC). (46) 

We will use / 33 « 1/20 with the understanding that it is uncertain by a factor of 2-3 either 
way. Note that / 33 is considerably larger, by about a factor of 200-700, than the value 
estimated in Ref. [T2|]. This results partly from lowering of m Vr from a few eV (used in 
Ref. fli~2| ) to about 1/20 eV (the SuperK value) and in part from interplay between the 
mixings in the Dirac and the Majorana mass matrices via the seesaw mechanism. The latter 
has the net effect of enhancing Mr w Bmjj/m u:j , for a given m U3 , precisely by a factor of 
B w 5, compared to what it would be without mixing. (Compare Mr w 10 15 GeV with its 
counterpart in Sec. II where, by ignoring mixing, we got m T R = (1 — 3) x 10 14 GeV.) The 
two effects together, - i.e., lowering of m U3 and increase through B - significantly enhance 
/ 33 = [Bmfj/m^/vR. 

It is interesting to note that the larger value of Mr arising because of mixing has a 
further implication. Using Mr w 10 15 GeV (corresponding to £ ~ 1), Vr — 2 x 10 16 GeV 
(for kr m 1) and / 33 w 1, one obtains: M = f 33 v R /M R k4x 10 17 GeV, which is the 
perturbative string scale [E0|. 



'Barring near-complete accidental cancellation between diagonal and off-diagonal contributions 
to m v . 
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Note that the net value of B (see Eq. (45)) depends on the parameters of the Majorana 
mass matrix as well as on a and e from the Dirac mass matrix of the neutrinos, which 
in turn are determined within 5*0(10) by the masses and mixings of quarks and charged 
leptons. That is why our expectations for proton decay are significantly affected by our 
understanding of the masses and mixings of quarks and charged leptons. 

The full set of contributions to the proton decay amplitudes (from all the operators in Eq. 
(41) involving all possible combinations of family indices) was obtained numerically using 
Mathematica and is listed in Appendix B.l. Calculations of a few representative (dominant) 
contributions to the amplitudes are exhibited in detail in Appendix B.2, where estimates 
of the amplitudes allowing for uncertainties in the relative phases of different contributions 
are presented. A general discussion of proton decay rate is given in Appendix C. Based on 
the result of these appendices, we now discuss the general constraint on the proton decay 
amplitude and thereby on the mass scale M e g- and M 16 tan 7 corresponding to the numerical 
estimate of the amplitude given in Appendix B. These constraints arise from existing lower 
limits on the proton lifetime. 

C. Constraints on proton decay amplitudes from the proton lifetime 

In Appendix C we show that with a certain (apparently) reasonable choice of supersym- 
metric spectrum, the lifetime of the proton decaying into neutrinos is: 

r _1 (p -> V T K + ) « (2.2 x 10 31 ) yrs x 

Here A(v) = A(u)/(2f), where A{v) is the strength of the four fermion proton decay am- 
plitude and / is the average wino-dressing function (see Eq. (93) in Appendix B and Eq. 
(111)-(114) in Appendix C). The quantity A(v) is simply the product of all the vertex factors 
in the wino-dressed Higgsino exchange diagram divided by the effective mass of the relevant 
color triplet Higgsino. For normalization purpose we can define A(v) slJ ^ = (X c X s 9q)/Mh c , 
where A, stand for the Yukawa couplings of the quarks at M\j. (For clarity of discussions, 
only the second generation contribution is kept here.) The quantity A(u) is the full ampli- 
tude, including the loop factor associated with the wino dressing. If we substitute A(z7) su ( 5 ) 
defined above into Eq. (47), we will reproduce the results given in Ref. |[34|| . (We have 
allowed for a factor of 4 enhancement in the lifetime relative to |j3~4"| , corresponding to an 
apparent slip by a factor of | in going from Eq. (3.7) to Eq. (3.8) of that paper.) 

Note that in writing Eq. (47), the short distance renormalization As of the d = 5 
operator in going from Mjj to Msusy (see Ref. ||34|| ) as well as the running factor to go from 
-^susy to 1 GeV have been included. As has a central value of about 0.67, which we shall 
adopt even for the 50(10) model. 

We have discussed in the appendix that typically proton would decay dominantly into 
V T K + and V p K + modes.0 In some cases the former supercedes the latter by factors of 2-5 



H + K° is also a likely prominent mode of the proton decay (see Sec. VI H) with a typical 



0.006 GeV 3 



(1/6) 



11 . 



1 2 



1 TeV 



2 x 10" 24 GeV" 
W) 



• (47) 
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in the rate, while the converse is true in some other cases (see Table. 1, Appendix C). We 
see from Eq. (47) that with A(j7 T ) k 2 x 10~ 24 GeV -1 , a reasonable "central value" for 
the partial lifetime r~ 4 (]9 — > V^K + )~ 1 is 2.2 x 10 31 yrs., which corresponds to the hadronic 
matrix element (3h = 0.006 GeV 3 (this is the central value quoted in lattice calculations 
||41||), (m^r/rriq) ~ 1/6 and rriq « 1 TeV. Now, adding the V li K + mode with a branching 
ratio R (-Raverage ~ 0.3, say), this is however 25(1 + R) times smaller than the empirical 
lower limit ||2| 

r -i(p -> 77^+) cxpt > 5.6 x 10 32 yr. (48) 

Thus, if the parameters have nearly their central values, the amplitude for the dominant 
of the two modes must satisfy the bound: AiV t K + ) < 2 x 10~ 24 GeV -1 /Vl + R ~ 4 x 
10~ 24 'GeV -1 / vT+~R, where £ = /x or r. Allowing that both /3# and the ratio of superpartner 
masses (m^/rriq) might well be smaller by factor of 2 (say) than the value quoted above, 
and that rriq could be (say) 1.4 TeV rather than 1 TeV, the theoretical value of the lifetime 
could plausibly increase by a factor of 32 compared to the "central value" 2.2 x 10 31 yr. This 
saturation of all uncertainties in the parameters, all in the direction so as to extend proton 
lifetime, strains credulity. If, nevertheless, one allows for such a variation in the parameters, 
the bound on the amplitude mentioned above will be relaxed by a factor of 5 to 6. Thus the 
empirical limit on proton lifetime leads to the rather conservative bound: 

A(pi) < 2.3 x 1(T 24 GeV^/VTTR . (49) 

This limit should be satisfied for every neutrino flavor V{ (assuming V fl K + and V T K + are 
comparable). For squark masses not exceeding about 1.5 TeV, we take Eq. (49) as a strict 
upper bound. It should however be noted that the lifetime depends quartically on the squark 
mass, so increasing rriq by a factor of 2 to about 3 TeV, holding fixed, would lengthen 
proton lifetime by a factor of 16. If true, this would relax the bound on A{p) quoted in 
Eq. (49) by as much as a factor of 4 to A(u) < 9 x 1(T 24 GeV~ 4 /y/i + R. Such a heavy 
spectrum (rriq ~ 3 TeV) for all three generations of squarks would however require severe 
fine-tuning of parameters in order to keep the vacuum expectation value of the light Higgs 
field at the electroweak scale. We shall assume a lighter squark spectrum (rriq < 1.5 TeV), 
for which there is no need for such an adjustment of parameters. In this case, the bound in 
Eq. (49) will have to be satisfied. 

D. Constraints on M e g and proton decay via standard operators 

In Appendix B, Eqs. (101) and (102) (or Table 1), we show that within the concrete 
5O(10) model, the dominant V T K + or V fl K + decay amplitude from the standard d = 5 
operator for cases I and II are given by: 

A(vK + ) std ~ 



branching ratio ~ 20-30%. This will however not alter our discussion here appreciably. 



2/4/M) 




'2.8 x 1(T 5 " 


M eff 




.1.2 x 1(T 5 . 



1 3 



(50) 
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There is an analogous expression for the new neutrino mass-related d = 5 operator, that will 
be discussed in the next subsection (E). The upper and lower entries in Eq. (50) correspond 
to cases I and II respectively. 

We now compare the upper limit (Eq. (49)) on the amplitude from proton decay searches 
against theoretical expectations based on the concrete 5*0(10) model. This leads to con- 
straints on M e g from the standard d = 5 operator (and on Miq tan 7 from the new operator). 
The upper bound (Eq. (49)) on A(Vi) applies to the net amplitude, which is given by the 
sum of the contributions from the standard (Eq. (50)) and the new operators (Eq. (51) 
below). For the sake of clarity, we will derive constraints on M c g (M 16 tan7) under the 
assumption that the standard (respectively, the new) operator dominates. Indeed, near- 
complete accidental cancellation between the two contributions is unlikely to occur for both 



v T K + and v^K^ modes. 



Using the net contribution from the standard operators to the amplitude given by Eq. 
(50), and the definition of A(Vi), we obtain (putting h 33 ~ 1/2) 



A{v T ) st d <*■ 



x 10~ 6 (l/2 to 3/2) . (51) 



Comparing with the empirical upper limit on A(z7j) (Eq. (49)) obtained as above, we get: 



M eff > 



4 

1.7 



x 10 18 GeV (1/2 to 3/2) . (52) 



Thus we see that M e $ has to be rather large compared to the MSSM unification scale of 
2 x 10 16 GeV in order that the standard operators may not run into conflict with the observed 
limits on proton lifetime. In effect, this reflects a net enhancement - by almost two orders 
of magnitude - of the standard d = 5 proton decay operators for realistic 5*0(10), compared 
to those in minimal SU(5), with low tan/? < 3 (see discussion in Appendix C). In the latter 
case, one need only require that the color triplet mass exceed about (2 — 3) x 10 16 GeV. 

Now, as mentioned in Appendix C, there are theoretically attractive mechanisms whereby 
the mass of 10'^, denoted by M w (see Eq. (38)), can be suppressed relative to the unification 
scale M\j. In this case, M c g = (Xa) 2 /Mio can be larger than Xa ~ M\j. Very large values 
of M cS ^> M\j could however lead to large positive corrections to a 3 (M z ), just from the 
doublet-triplet mechanism, above and beyond the value expected on the basis of simple 
coupling unification. For the doublet-triplet splitting mechanism described by Eq. (39) the 
shift in a^{mz) from this sector alone is found to be (see Appendix D) 

a / \ 1 M™z)] 2 9 /M eff cos7\ 

Aa 3 (mz)\ DT = {-^~ ) ■ ( 53 ) 

This generalizes the expression given in Ref. fll~2|j43|| , where the MSSM Higgs doublets were 



assumed to be contained entirely in 10 H and not in 16#, corresponding to cos 7 = 1. The 
argument of the logarithm in Eq. (53) is simply the ratio: (product of the three color triplet 
masses)/ (product of the two superheavy doublet masses xMu). From the determinant of 
Eq. (39) we see that the product of the color triplet masses is equal to M e sM w MiQ, while 
the two heavy doublets have masses given by M 10 and M 16 / cos 7. Note that the second 
heavy doublet has a mass larger than M 16 . This is the reason for the presence of the cos 7 
factor in Eq. (53). 
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To evaluate the RHS of Eq. (53), we use the lower bound of M e fj that is suggested by 
proton lifetime constraints, viz., M eff > (2-6) X 10 18 GeV or (0.8-2.5) x 10 18 GeV (see Eq. 
(52)), and the MSSM unification scale ofMu~2x 10 16 GeV. We should also specify the 
value of cos 7. It is obtained in terms of tan f3 as follows. From m t ~ h 33 (10n) n = h 33 v u and 
m b ~ h 33 (10 H ) d = h^cos^Vd, we have m t /m b ~ {v u /v d )(l/ cos 7). Inserting m t /m b ~ 60, 
we thus obtain cos 7 ~ (tan/3/60). (This can also be expressed as tan (3 tan 7 ~ m t /m b , 
which is valid for tan 7 > 3.) We see that cos 7 is a small number (« 1/30 — 1/6) for small 
and moderate values of tan f3 (~ 2 — 10). 

Now, let us recall that, in the absence of unification-scale threshold and Planck-scale 
effects, the MSSM value of a 3 {m z ) in the MS scheme, obtained by assuming gauge coupling 
unification, is given by a 3 {m z )\MSSM = 0.125 — 0.13 0,0. This is about 5-10% higher 
than the observed value: ot 3 {m z ) = 0.118 ± 0.003 |26|]. Substituting the central value, 
a 3 (M z ) = 0.118, and the MSSM unification scale of M n = 2 x 10 16 GeV, one obtains for 
cos 7 = 1/20: Aa 3 (m z )\ DT ~ 0.0046 - 0.0077, for M cS « (2 - 6) x 10 18 GeV. Thus the 
constraint on M e s from proton lifetime for case I amounts to having in MSSM a net value 



« 3 K)U t = a 3 (M z )^\ MSSM +Aa 3 (M z )^ M +A' 3 ~ (0.132-0.135)+A 3 , for cos 7 = 1/20, 
where A 3 denotes other unification scale threshold and Planck scale effects evaluated at the 
electroweak scale. Including Aa 3 (m z )rjT and A 3 , the net theoretical value of a 3 {m z ) is 
given by a 3 (m z ) nct = ct 3 (m z )^ SSM + Aa 3 (m z ) DT + A' 3 . Since a 3 (m z )^ SSM is higher than the 
observed value, and Aa 3 (m z ) OT is relatively large (depending on M ef j 00S7) and is positive 
(see below), one would need a net appropriately large negative value of A 3 so that a 3 (m z ) ne t 
may agree with the observed value. By varying the parameter T = M eff cos 7/ (2 x 10 18 GeV) 
appropriately so as to include the range M e ff > (1 — 6) x 10 18 GeV, obtained from proton 
lifetime constraint, we can evaluate Aa 3 (m z ) net and thereby assess the needed value of 
A 3 = a 3 (m z ) ohs - a 3 (m z ), where a 3 {m z ) = a 3 (m z ) M ssM + Aa 3 (m z ) DT . This is shown in 
the Table below. 



T 


1/60 


1/40 


1/30 


1/20 


1/10 


1/3 


Aa 3 (m z ) DT 
a 3 (m z ) 


0.0014 
0.1284 
-8.8% 


0.0026 
0.1296 
-10% 


0.0034 
0.1304 
-10.5% 


0.0046 
0.1316 
-11.5% 


0.0066 
0.1336 
-13.2% 


0.0100 
0.1370 
-16.1% 



Here 5' 3 = A' 3 /a 3 (m z ) ohs = (a 3 (m z ) obs - a 3 (m z ))/a 3 (m z ) ohs . In above, we have used a 
reasonable lower limit on ct 3 (m2)^[ gsM = 0.127 for nig < 1 TeV [42] and have used the 
central value for a 3 (mz) b s = 0.118. Allowing for a 3 (m z ) ohs = 0.118 ± 0.003 would amount 
to adding nearly ±2% change to 5' 3 . For concreteness, we will quote results for central value 
of a 3 (mz) k> s , but bear in mind the possibility of the ±2% change in 5' 3 . Note that the 
variation of the parameter T used above include proton lifetime constraints (Eq. (48)) of 
M efT > (2 - 6) x 10 18 GeV for case I, with cos 7 w 1/60 to 1/9, and M eS > (0.8 - 2.4) x 10 18 
GeV for case II, with cos 7 ~ 1/60 to 1/4. 

The Table shows that if coupling unification should hold, one must assume, for the case 
of MSSM embedded in 5*0(10), that other unification scale threshold and Planck scale ef- 
fects (excluding Aa 3 (m z )rjT), denoted by A 3 , provide a substantial negative contribution 
to a 3 (m z ), typically varying between —8% to as much as —16%, as evaluated at the elec- 
troweak scale. Before discussing the feasibility of such large negative correction, it should 
be emphasized that the presence of cos 7 inside the logarithm of Eq. (53) has played a 
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significant role in diminishing the positive threshold correction to a 3 (mz)- In its absence, 
with M c fj > (1 — 2) x 10 18 GeV, negative threshold corrections exceeding even (18-20)% 
would have been required from the other unification scale threshold and Planck scale effects. 

We now show that it is reasonable to expect the net threshold correction from other 
effects, denoted by A' 3 , evaluated at the electroweak scale, to be negative in sign, and to 
have magnitude no more than about 6 to 10%, if one confines oneself to low dimensional 
Higgs multiplets such as 45,16,16 and 10 (as we do). First, let us note that since the 
coupling and the correction get enhanced by nearly the same factor (« 3 for a 3 ) as they run 
from Mu to mz, a -10% correction to a 3 at m z roughly corresponds to a -3.3% correction 
at the unification scale. The precise value of the correction depends on the details of the 
model, including the nature of the Higgs system that causes mass splittings within complete 
SU(5) multiplets. Fortunately, as previously observed by several authors (see e.g. p7| , 
HH), for Higgs multiplets that are not too large (such as 45 and 16 of 5*0(10) as in our 
model) and with the split masses (M a ) of the sub-multiplets belonging to a given 5*0(10) 
multiplet being within a factor of < 2 — 10 of Mu either way, the net threshold correction to 
a 3 (mu), evaluated at the unification scale[] is typically no more than about 1 to 2%, and the 
corrections from a given sub-multiplet can be either positive or negative depending upon 
whether M a is greater or smaller than Mu, and also on the SU(3) x SU(2) x U(l) quantum 
numbers of the sub-multiplet. 

To get a feel for the magnitude of the correction, consider (for illustration only) the case 
of minimal supersymmetric SU(5) (with 24, 5, 5 of Higgs). In this case, owing to the color 
triplets (H 3 , -£%) in (5, 5), which are separated from the doublets (albeit by fine-tuning) and 
acquire unification scale masses, a 3 (mu) receives a correction Aa 3 (Mu) H3 = a 3 (Mu)(e 3 )# 3 
at the unification scale, where [42] (e 3 )# 3 = [3a un ifln(M# 3 /Mu)]/57r m +(0.5 — 1.7)% for 
M H3 /Mu ^2 to 10, and a uni f ~ 0.04. Note this correction is expected to be positive because 
proton lifetime constraint suggests M# 3 > Mu, and even for the extreme value of the mass 
ratio of 10, it is less than 2%.Q 

In the context of the minimal Higgs system (like ours) that utilizes the VEVs of just the 
45, 16 and 16 of Higgs to break 5O(10) to the standard model gauge symmetry, we can 
identify two (rather definite) sources of negative contribution to a 3 {mz). These arise from 



9 Note that the threshold corrections owing to doublet-triplet splitting in SO (10) (discussed in 
Sec. VI. A) is a special case, because the combination M c g cos 7 = (Aa) 2 COS7/M10, that enters into 
the logarithm in Eq. (53), does not represent the masses of the color triplets (which is of order 
Aa). This arises because of the cos 7 factor in the logarithm, and also because Miq representing 
the mass of the complete multiplet of 10' (see VI. A) can be much smaller than M\j (as needed 
from proton decay constraint, Eq. (52)). This is the reason that for M e g cos 7 > 2 x 10 18 GeV, 
Aa3(mz)DT can lead to large positive corrections to o 3 at mz- 

10 One can obtain large negative contribution to a 3 (mz) in SU(5) by introducing large Higgs 
multiplets, such as 75 and 50 [42] which are used in the missing partner mechanism for doublet- 
triplet splitting in SU(5). This however is a special case associated with large multiplets and does 
not apply to lower dimensional multiplets like 45, 16, 10 of SO(10). Large multiplets such as 75 
of SU(5) do not "explain" the observed unification of couplings, but rather accommodate it. 



26 



mass-splittings within (a) the gauge multiplets, and (b) the Higgs sub-mult iplets 45h- 

First consider the gauge multiplets. Representing the VEVs of 16# and 45h by c and a 
respectively (see Appendix D), one finds: 

Aa 3 (m z ) gaugc ~ (+0.002, -0.00026, -0.0015, -0.0028, -0.0031, -0.0030, -0.0018, -0.001) 
for 

p= (0.1,0.2,0.3,0.5,0.7,1,2,3), 

where p = 2c/ a. Thus the contribution can be positive for sufficiently small p, but for most 
of its range, p > 0.2, the contribution is negative, with (5 3 ) ga uge = Aa3(m z ) gauge /a;3(m z ) bs 
varying from —1.3 to about —2.6%, for p varying from 0.3 to 2.0. Note that the maximum 
of —2.6% evaluated at m z corresponds to a correction of about —0.8% at the unification 
scale, in accord with expectations. 

The contribution to 0:3 (mz) from the splitting of 45h has been evaluated in Appendix 
D (see Eq. (119)). As shown there, if one ignores the coupling of 45h to other multiplets, 
the allowed superpotential has a simple form: W = Mi(45h) 2 + k(45h) 2 /M. One can then 
argue that the masses of the submultiplets are characterized by Mi « n(My/M) ~ 10 _2 Mu, 
if k ~ 1 and M ~ M Planck . One then obtains Aa 3 (m z ) 4 5 H ~ —0.0045. 

Thus, at least in the simplest approximation, the contribution from 45h to (mz) is 
negative and is about -4% at the electroweak scale. Of course, the above superpotential 
is a simplification because one needs to couple 45h to other fields that acquire unification 
scale VEVs, like 16h, 16h, together with possible additional vector-like pairs (16v, 16v) 
without VEVs to avoid uneaten pseudo-Goldstone bosons. But such couplings do not alter 
the basic feature that Mi -C Mu (see [B7[|). 

Thus we see that with the minimal Higgs system (i.e., with just 45h, 16h and 16h 
acquiring VEVs), there are good reasons for the unification scale threshold correction to 
otz(mz) from our "other effects" denoted by A 3 , to receive negative contributions from the 
gauge sector of —1.3 to about —2.6% (depending on p), and plausibly about —3 to —4% from 
the Higgs multiplet 45h at the electroweak scale. These two together can easily combine to 
yield about —4.5 to —6.6% correction to a 3 (m z ). 

Having identified two rather definite sources of negative threshold corrections to a^(m z ), 
the remaining other contributions (excluding Aa^(mz)DT) arise through splittings within 
16h, 16h, 16v and 16v (which would be induced through < 45h >)• Contributions to 
each of these multiplets can be either positive or negative depending on the parameters in 
W. Since the induced splittings within these multiplets will be relatively small, as also their 
sizes (for mass ratios < 2), the magnitudes of these corrections from each of these multiplets 
to a 3 (niz) can be estimated to be less than about 1%. Thus the combined correction from 
16h, 16h, 16v, 16v is expected to be less than about 2 to 4% (positive or negative) at the 
electroweak scale. 

Finally, it is worth noting that with only 45h, 16h and 16h acquiring VEVs (as in 
our model), contribution to a^(m z ) from Planck scale physics through effective operators 
F^Fn^bu/M will vanish because of antisymmetry in the 5*0(10) contraction. 

Thus, at least for the minimal Higgs system utilizing the VEVs of only 45h, 16h, 16h, 
we see that the net "other" threshold correction to 0:3 (mz) from unification and Planck 
scale effects, excluding Aa 3 (m z )r)T but including those from (a) the gauge sector (—1 to 
-2.6%), (b) from 45 H (-3 to -4%), (c) collectively from 16 H ,16 H , 16 v ,16 v (±2 to 4%), 
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and (d) from Planck scale effects (« 0%), is very likely negative. But after adding all these, 
it seems plausible to assume that the magnitude of this net other contribution is not more 
than about 8 to 10%, at the electroweak scale. 

This assumption in turn implies (see the Table above) that T = M c r cos 7/(2 x 10 18 GeV) 
is bounded from above by about 1/30, or conservatively by 1/20. If we add a +1.5% 
correction to the entries for 5' 3 in the Table to allow for otz{mz) = 0.118 + 0.0017, then 5' 3 
(reduced as above) is < 10%. This in turn leads to the upper bound: 

For a given cos 7 ~ tan/3/60, this upper bound yields an upper limit on M e g and thereby 
an upper limit on proton lifetime. In particular, if we assume tan/5 > 2, and thus 
cos 7 > 1/30 we obtain M eff < 3 x 10 18 GeV. Putting tan/3 = (2,3,6,8,10) and 
thus COS7 = (1/30,1/20,1/10,1/7.5,1/6), the upper bound given above yields (M e g/10 18 
GeV) < (3,2,1,0.75,0.6). Combining this upper limit on M eff (which holds for both 
cases I and II) arising from the threshold corrections to a^iviiz) with the lower limits on 
(M efr /10 18 GeV) > (2-6) for case I and > (0.8 - 2.4) for case II, given by the proton 
lifetime constraint (Eq. (52)), we see that the two limits would be in conflict with each 
other if 00S7 > (1/20, 1/7.5) - i.e., if tan/3 > (3,8) - for cases (I, II). These considerations 
based on Aa^m^) and the experimental limit on proton lifetime suggest that rather small 
values of tan/3 -i.e. tan/3 < 3 (case I) and tan/3 < 8 (case II) are favored (but see below). In 
either case, for the MSSM, assuming tan/3 > 2, by demanding accurate coupling unification 
we obtain an upper limit on M e g given by 

Mcfr < 3 x 10 18 GeV (55) 

On the other hand, the proton lifetime constraint (Eq. (52)) implies that M e g- must 
exceed (2 to 6)xl0 18 GeV for case I and (0.85 to 2.5) xlO 18 GeV for case II, where the range 
in M e fj corresponds to the uncertainty factor B =(1/2 to 3/2) in the amplitude (see Eq. 
(51) and Appendix B) in a correlated manner. For instance, for M e g < 3 x 10 18 GeV (rather 
than 6 x 10 18 GeV) B can vary only between (1/2 to 3/4) for case I, while for case II, B can 
still vary between (1/2 to 3/2). Using this range and M e g < 3 x 10 18 GeV, we can obtain a 
lower limit for the proton decay amplitude (given by Eq. (51)): 

Mv K+) , > f (T X lO-^GeV^Xl/e to 1/4) 

A{u T K ) std > [ (3 x 1(r24 Ge y-i )(1/6 tQ 1/2) J (56) 

Substituting into Eq. (47) and adding the contribution form the second competing mode 
V„K + with a typical branching ratio R ~ 0.3, we obtain 



T-\vK + ) std < 



(3 x 10 31 yrs.)(1.6 to 0.7) 
(6.8 x 10 31 yrs.)(4 to 0.44) 



(32 to 1/32) (57) 



If tan/3 < 1.5, the top Yukawa coupling will blow up before reaching [i = My. Within the 
standard supergravity spectrum of super symmetry breaking, there are indications that tan (3 > 2 
from LEP constraints. 
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Here the uncertainty (32 to 1/32) corresponds to the uncertainty in /3#, ( m wl m <i) anc ^ m <?> 
by factors of 2,2, and \/2 respectively, either way, around the "central" values reflected in 
Eq. (47). Thus we find that for MSSM the inverse partial proton decay rate should satisfy: 



T-\p^vK + ) stA < 



3 x 10 31±L7 yrs. 



6.8 x 10 31 -i'-s yrs 

(MSSM) (58) 



< 



1.5 x 10 33 yrs 



7 x lO^yrs 



The upper limit in Eq. (58) essentially reflects the upper limit on M e g, while the remaining 
uncertainties of matrix elements and spectrum are reflected in the exponents. We stress 
that the upper limit on lifetime exhibited in Eq. (58) has arisen by first taking an upper 
limit of 10% on the net unification scale threshold correction to — a% (excluding the contri- 
bution from the doublet-triplet splitting), which appears to us generous. Furthermore, the 
amplitude limit exhibited in Eq. (56) is obtained only if the uncertainties in the amplitude, 
Ph, ( m w/ m q) an d rriq all go in the same direction to about their extreme values to extend 
proton lifetime, so it is a generous limit. 

Before discussing the contributions of the new operator to proton decay, we wish to note 
an interesting possibility, that can be relevant especially if tan/? is large. For large values 
of tan/? « 20 - 30, corresponding to cos 7 « 1/3 to 1/2, with M eS > (2 - 6) X 10 18 GeV 
(which is needed for case I to satisfy the proton lifetime constraint), one would have T = 
M efr cos7/(2 x 10 18 GeV) > 1/3 to 3/2. In this case, in the context of MSSM, a negative 
threshold correction (5' 3 ) exceeding 16 to 20% in magnitude would be required from other 
sources (see the Table given above for 5' 3 as a function of T). As we said before, although 
such a large negative correction, together with a matching positive one, is in principle 
possible (see e.g. Ref. [43]), it diminishes the luster of the observed agreement of simple 
coupling unification, by making it appear somewhat fortuitous. In this connection it is 
noteworthy that extra vector-like matter - specifically a 16 + 16 as proposed in the so-called 
ESSM (Extended Supersymmetric Standard Model) [[4(J - at the TeV scale could greatly 
ease this problem, and thereby allow large values of tan/3, while leaving our discussion of 
ordinary fermion masses essentially unaltered. In this case, a un if is raised to nearly 0.25 
to 0.3, compared to 0.04 in the MSSM. Owing to increased two-loop effects, the scale of 
unification M\j is raised to (1 — 2) x 10 17 GeV P6]J47[], while a 3 (Mz)°\essm is lowered to about 



0.112 — 0.118. With increased My the correction Aa 3 (M z )\ DT is also lowered. As a result, 
even for tan/3 ~ 20 — 30, i.e., COS7 ~ (1/3 — 1/2), and M cS »6x 10 18 GeV, one obtains 
for the case of ESSM, the net value of a 3 (m z ) = ct 3 (mz)°\ESSM + Ac^m^^r + A3 k 
(0.123 - 0.124) + A3. Thus for ESSM embedded in 5O(10), unification scale threshold 
corrections from "other" sources to a 3 (m z ) (represented by A 3 /a 3 (m^)), though negative, 
need be no more than 5% in magnitude, even in the rather extreme case tan/3 = 30, with 
M e g = 6 x 10 18 GeV. As shown above, for the minimal Higgs system a net negative threshold 
correction to a 3 {mz) of this magnitude is not unexpected. As an added feature, we note that 
since in the ESSM case M e g can be rather large, ~ (4—8) x 10 18 GeV, without requiring large 
negative threshold corrections from other sources, we can be compatible with the observed 
limit on proton lifetime for the central values of uncertainties in the parameters describing 
the spectrum of the supersymmetric particles and the relevant matrix elements. 

To be specific, for M c r « 6 x 10 18 GeV and tan/3 = 30, which corresponds to — 5' 3 = 4 
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to 5%), we obtain (compare with Eq. (58)): 
T- 1 (uK + ) std 

where the last factor arises from allowing (m^/m^mg, and /3# to vary by factor of 2,2, 
respectively about their "central" values. The upper and lower entries correspond to cases 
I and II respectively. Note that allowing for a modest factor of 1 to 4 jointly from the two 
brackets for case I (and 1 to 2 for case II), which correspond to nearly central values of 
the parameters mentioned above, keeps the theoretical value of proton lifetime within the 
experimental limit. On the other hand, allowing for a factor of 20 jointly from the two 
brackets (for either case I or case II), the proton lifetime in the ESSM case typically lies in 
the range (1 — 5) x 10 33 yrs, which should be accessible. 



1.2 x 10 oz yrs. (1.6 to 0.7) 
2.7 x 10 32 yrs. (4 to 0.44) 



(32 to 1/32), ESSM 



(59) 



E. Constraint on M±e tan 7 and proton decay via the new operator 

The decay amplitude for the new operator for the leading mode (which in this case is 



v iJL K + ) is given by 
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y Mi6 tan 7 
e Q/ 3 7 (GPw /3 )(s 7 z/ 3 ) 



(3x 10- 6 )(l/2to2) [f(t,d) + f(t,l)]x 



(60) 



as shown in Appendix C. This is true (approximately) for both cases I and II. 
From the definition of A(v) we obtain: 



/ijncw 
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Mi6 tan 7 



(3 x 10~ 6 )(l/2 to 2) 



(61) 



The first bracket P = [/ 33 /i 33 /(Mi 6 tan7)] may be evaluated by using results of Sec. VI A 
and B, as follows. Putting j 33 = f 33 VR/M (see Eq. (46)), Mi 6 tan7 = X'vr (see Eq. (40), 
and /i 33 ~ 1/2, we obtain: 



P = (/ 33 /M)(l/2A') . 



(62) 



Here M stands for the mass scale that characterizes the strength of the effective non- 
renormalizable operators (see Sec. III). Using Mr = f 33 v R /M and the result that Mr ^ 
Biml/m^) ~ 5[m 2 (Mu)/(l/20 eV Q ^ 10 15 GeV/C (see Eqs. (44) and (45)), we find:Q 



/33/M = Mr/ v r « (4Cr[4 x 10 17 GeV] 



1-1 



(63) 



12 We note that the 5O(10) contraction for the coupling contributing to neutrino masses is not 
the only one that contributes to proton decay. We assume that the two possible contractions are 
comparable in strengh. 
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Here kr and ( denote uncertainties in vr and m Vz - i.e., vr = kr(2 x 10 16 GeV) and 
m U3 = ((1/20) eV. A priori, we expect vr ~ My and thus /tR ~ (1/2 — 3). Since < 16# > 
breaks S'O(IO) to SU(5), vr =< 16# > may in fact be somewhat larger than My - that 
is, k r > I. From the SuperK results, we have ( 1.5 — 1/1.5. In Sec. Ill, we have 
discussed the appropriateness of the characteristic mass M being either the Planck or the 
string scale, subject to the presumption that the effective coupling / 33 for the third family 
is nearly maximal (~ 1). Given some uncertainty in this regard, a good choice appears 
to be M = 10 18 GeV, which is intermediate between M P i anck and M string . We expect that 
such an intermediate value of M should represent fairly well either choice M = M Planck or 
M = M string , given that we allow a range in k r ( and thereby in / 33 . Setting M = 10 18 GeV, 
Eq. (63) yields: / 33 ~ 2.5/(k 2 r (). Noting that values of / 33 1 are implausible, we thus 
expect k 2 r ( > 1 in accord with the remark mentioned above (rather than < 1). Under the 
presumption that / 33 is maximal, it therefore seems quite reasonable to assume that k 2 r ( ~ 
(1 to 5) 2.5(1/2 to 2), which corresponds to / 33 ^ (2 to 1/2). Substituting this range for 
k r ( into Eq. (63) and (62), we obtain 

P « (5 x 1(T 19 GeV _1 )(l/2 to 2)/A' 

. Now, on the one hand, the validity of the perturbative treatment requires A' < 1. One the 
other hand, given that M 16 tan 7 = X'vr and that from considerations of standard operators 
for proton decay, we have the constraint that tan (3 < (3, 8) and thus tan 7 > (20, 7) for cases 
(I, II) in MSSM, the choice of A' < 1 will tend to make M 16 <C vr, which is implausible. We 
therefore take A'wl. Substituting this into P and in turn into Eq. (60), we obtain 

A(V^K + ) 11CW « (1.5 x 10~ 24 GeV _1 )(l/4 to 1.3) (64) 

Here the upper end of the uncertainty has been restricted to conform to the limit on proton 
lifetime, see Eq. (49). 

Comparing with Eq. (56) we see that the contributions of the new and the standard 
operators (for M e g (2 — 3) x 10 18 GeV) to the proton decay amplitude are comparable to 
one another. Since there is no reason to expect near cancellation between them (especially 
for both V^K + and VK + modes), we expect the net amplitude (standard + new) to be in the 
range exhibited for either one. Assuming that the new operator dominates and substituting 
Eq. (64) into Eq. (47), we obtain: 

r- 1 (i7 J ftT + ) new ss (3 x 10 31 yrs)[16 to 1/1.7]{32 to 1/32} (65) 

In this estimate we have included the contribution of the V T K + mode with a typical branch- 
ing ratio R ~ 0.4 (see Appendix C). Here the second factor, inside the square bracket, 
reflects the uncertainties in the amplitude (see Eq. (64)), while the last factor corresponds 
to varying fin, {^w/ 171 ^) anc ^ m <i aroun d the central values reflected in Eq. (47). 

With a net factor of even 20 to 100 arising jointly from the square and the curly brackets, 
i.e. without going to extreme ends of all parameters, the new operators related to neutrino 
masses lead by themselves to proton decay lifetimes ~ (0.6 — 3) x 10 33 yrs. Thus if the new 
operators were the only source of proton decay - i.e., if the standard d = 5 operators were 
somehow absent - one could be comfortably compatible with existing limits, but optimistic 
regarding future observation. We now briefly elaborate on this possibility. 
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F. The neutrino mass related operator as the sole source of proton decay 



As we have seen, straightforward minimal embedding of the MSSM in 50(10), with 
informed hypotheses about the fermion-Higgs couplings, leads to contributions from the 
standard d — 5 operators are disturbingly large, if the relevant parameters have nearly their 
central values and/or if M c g < 2 x 10 18 GeV. This is especially true for case I in the MSSM 
(see Eq. (57)- (58)). One is therefore motivated to wonder whether the standard operators 
might not be present. This possibility is realized, if the higher gauge symmetry is G224 (or 
C2113 = SU{2) L x I 3R x (B — L) x SU(3) C ) rather than SO(10). Such gauge symmetries have 
been proposed to appear in solutions of string theory (see Ref. for G224 and Ref. for 



G2113). Such possibilities retain some attractive features of SO(10), notably the unification 
of quark-lepton families into single multiplets including the right-handed neutrino v R as 
indicated for the neutrino seesaw. One sacrifices a simple group-theoretic explanation for 
the observed unification of coupling, but if the models derive from an underlying string 
theory, one might still expect such unification at the string scale. Plausible mechanisms 
to reconcile the string and the MSSM unification scale have been proposed |5(|. The color 
triplets related to the electroweak doublets, which generate the standard d = 5 proton decay 
operators, need not exist. The standard source of d = 5 operators can be absent in such 
models pT|j . 

It is possible that after projecting out appropriate fields the couplings of those that 
remain reflect the original higher symmetry^. This is what is believed to occur for the gauge 
couplings, as previously mentioned. If it also occurs for the Higgs superpotential couplings, 
our considerations on fermion and neutrino masses in Sec. II-IV and their relationships to 
proton decay will remain valid. 

G. Baryon number violation from the RRRR operator 

So far we have focused on the charged wino dressing of the effective superpotential of 
Eq. (41). For small values of tan f3 and fi, this gives the dominant contribution to the proton 
decay amplitude. However, if /ztan/3 is large, (> 2 TeV), dressing of the effective baryon 
number violating operator involving only the right-handed fields, u c u c d c e c , by the charged 
Higgsino can become important H53|-[55|]. The ratio of this amplitude to the usual wino 
contribution scales as tan j3(fj,/m^r) (for rriq 3> m^). In minimal supersymmetric SU(5), the 
Higgsino dressing becomes more important than the wino dressing when tan/3 > 9(111^/(1) 
|55| . This estimate takes into account the differences in the renormalization of the (RRRR) 
compared to the (LLLL) operator owing to running from M\j to Msusy, and the difference 
between their strong matrix elements. We now show that, by contrast, in the S'O(IO) model 
the RRRR operator is strongly suppressed relative to the LLLL operator, and that they 
will be comparable only if (//tan ft/m^) ~ 500. 



13 For example, in a class of string solutions leading to Cr2ii3> the cubic level top and v T Yukawa 
couplings are claimed to be equal at the string scale despite SU (4)c breaking 
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In the 5*0(10) model, once the quark and lepton masses and mixings are fixed, the 
strength of the RRRR operator is determined by the corresponding effective superpotential 
(analogous to W^) which is given by: 

WS ) = -M;i[2{u c T HV'e c ){u c T HKd c ) + (u c T HV'e c )(u c T AKd c ) 

- 3(u c T AV'e c ){u c T HKd c ) - 3(u c T AV'e c ){u c T AKd c ) 

- tan 7 (w c T HV'e c )(u c T GKd c ) + 3{u c T AV'e c ){u c T GKd c )} 

+ M' 1 {u c T FV'e c ) {u c T GKd c ) . (66) 

Here 

(H,A,G,F) = V^(h,a,g,f)V u , (67) 

with V u the unitary matrix that rotates the u c fields u c ^ = VuU ^ from the gauge basis 
to the mass basis, V e is the unitary matrix that similarly rotates the right-handed electron 
field, and V' = V^V e . K is the right-handed analog of V C km, K = V^V d . 

The contribution to proton decay amplitude from Eq. (66) is estimated as follows. Since 
a charged Higgsino is involved in the dressing, the internal scalars have to be from the third 
generation (other diagrams will be suppressed by small Yukawa couplings). This uniquely 
picks out the Ir and f R as the internal scalars. The external quark fields are then fixed to 
be u and either a d or an s. This suggests that the combination of indices (ij)(kl) in Eq. 
(66) must be (33) (11), (33)(12), (13)(32) or (13) (31) . Among these four, the combination 
(33) (12) can proceed without utilizing any of the right handed mixing angles. (Though the 
right-handed mixing angles in the (23) sector of both u and d are ~ 0.2, not terribly small.) 
We find the amplitude, after Higgsino dressing, from this dominant contribution to be 

A[(u c s c Y(du T )] ~ M^hl^'Vu + (M 16 tan 7 )- 1 /i 3 3iW^ (68) 
w 6.6 x 10~7 M cfr + 1.3 x 10- 24 GeV _1 (l/2 to 2) . (69) 

Other contributions are not much bigger. In going from Eq. (68) to (69), we have used 
h 2 33 w 1/4, i ~ 4.4 x 10~ 3 (see Sec. IV), V td ~ 6 x 10~ 3 and P = f 33 h 33 /M 16 tan 7 ~ 
5 x 10~ 19 GeV _1 (l/2 — 2) (see the discussion following Eq. (63)). It is understood that the 
relative sign (phase) of the two terms in Eq. (68) and (69) is arbitrary. The full amplitude 
A[(u c s c Y (dv T )} will be obtained by multiplying the above expression by a loop function 
analogous to the function for the wino dressing (see Eq. (92)). There are two differences 
in this function however: (i) The factor (a 2 /47r) will be replaced by (A t A r )/(167r 2 ), where 
the Yukawa couplings are to be evaluated at the momentum scale Msusy- (h) The mass 
parameter is replaced by /i. The presence of A T brings in a tan/3 dependence in the 
Higgsino dressing relative to the wino dressing. (In minimal super symmetric SU (5), the wino 
dressing has a tan (3 dependence, but the Higgsino dressing will have a (tan/9) 2 dependence.) 

Comparing the first term of Eq. (69) with the standard LLLL amplitude (see Eq. 
(51)) and making the replacements as above, we obtain: ^4_RR_Ri?(std)/Aiiii(std) ~ 
[h t h rf i/(g 2 2 m w )}V ts x {6.6 x 10- 6 }/{[7,3] x 10- 6 }(l/2) » (// tan P/m^)V ts [l/100, 1/48], 
where we have put h t (m t ) xs 1 and h T (m t ) xs m T /v D xs tan/3/100. Substituting V ts ~ 1/25, 
we thus see that even for relatively large values of (fxtan/3/m^ 100, the RRRR amplitude 
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is smaller than the LLLL amplitude by a factor of 25 — 12. Thus, unlike the case of 5C/(5)Q, 
in the SO (10) model developed here, the RRRR operator can safely be ignored compared 
to the LLLL operator, as long as (/itan/^/m^) < 200 (say). 

The charged Higgsino dressing results in decay modes of the proton containing final state 
neutrinos, not charged leptons. Apart from the direct charged Higgsino exchange, there are 
diagrams involving ti — mixing and ft — Tr mixing (the latter being proportional to 
tan/3) followed by the exchange of charged wino. This contribution is smaller than the 
direct charged Higgsino exchange by roughly a factor v u /M S u SY - Similar arguments apply 
to mixed contributions involving charged Higgsino-wino mixing. 



H. Charged lepton decay mode 

In minimal supersymmetric SU(5) and many of its variants, charged lepton decay of 
the proton is suppressed as long as tan/5 < 20 or so. This is because of a GIM-type 
cancellation in the wino dressing diagrams, which brings in a suppression factor in this 
amplitude, proportional to the small w-quark mass. Gluino dressing can lead to charged 
lepton decay of the proton, but owing to flavor conservation of the primary gluino vertex 
and the flavor antisymmetry of the effective superpotential, this contribution is suppressed 
for small values of tan/? (< 20). For large values of tan/5 (> 20), flavor mixing in the 
up-squark sector becomes significant and the gluino graph begins to be important f)6| . 



However, for minimal supersymmetric SU(5), such large values of tan/5 (> 20) are highly 
disfavored owing to limits on proton lifetime (see for example Ref . ) . Even if we ignore 
this difficulty, for large tan/3 > 20 dressing of the RRRR operator by charged Higgsino 
becomes the dominant source of proton decay in SU(5) (if /V m w> > 1), and the neutrino 
modes dominate anyhow. So in minimal supersymmetric SU(5), and many of its variants, 
charged lepton decay mode of the proton is highly suppressed relative to the neutrino modes, 
for all values of tan /3. 

The situation is different in our S'O(IO) model, for two reasons. First, the contribution to 
p — > {i + K° arising from the standard d = 5 wino dressing diagram in this model is not small. 
Since we have a realistic spectrum of quark and lepton masses, especially with m s ^ m M and 
rrid 7^ Tn e at Mu, the wino contribution does not experience a GIM cancellation. Second, the 
new d = 5 operator related to neutrino masses lends comparable strength to the neutrino 
and the charged lepton modes. 

In evaluating the expected branching ratio for the fi + K° mode - i.e., B([i + K°) = 
T(jj, + K°)/[T(V T K + + T(V fl K + ) + r({i + K )} - let us concentrate on case I. Very similar 
results hold for case II. The amplitudes for p — > fi + K° arising from the standard and the 
new operators for case I are given by (see Appendix B, Eq. (107) and (109) or Table 1): 



The main reason why the standard RRRR operator is suppressed, relative to the LLLL operator, 
in SO(W) but not in SU(5), is simply that the standard LLLL operator is effectively enhanced 
by about two orders of magnitude in SO(W), for M e g = Mq- See the remarks following Eq. (114) 
in Appendix C. 
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h 2 

A{n+K°) atd » ( T p-)(3 x l(T 6 )(l/2 to 2) w (0.30)(l/2 to 2) x lO-^GeV" 1 (70) 

•Meff 

i(/i + K°) ncw » P(HT 6 )(l/3 to 2) » (0.5)(l/6 to 4) x KT^GeV -1 (71) 

In Eq. (70) we have inserted /ig 3 ~ 1/4 and an average value of M e g- m 2.5 x 10 18 GeV 
(see discussions following Eq. (55), showing M e s ~ (2 — 3) x 10 18 GeV for case I). In Eq. 
(71), we have inserted the value of P = / 33 W(Mi 6 tan 7) w (5 x Kr 19 GeV _1 )(l/2 - 2). 
obtained in Sec. VI E. The standard and the new amplitudes for the T7 T K + and the V IJL K + 
modes (for M eS ~ 2.5 x 10 18 GeV) for case I are given by (see Eq. (101), (102), (105) and 
(106) and Table 1 in Appendix B): 

^(^)^{ (1 (i!lor^ 2) }x«v- cb, 

Note that the upper ends of the ranges shown for the different amplitudes have been re- 
stricted to conform with the limit on proton lifetime (see Eq. (49)). Assuming that the 
relevant matrix elements for the p> + K° and VK + modes are comparable (but see below), it 
may be inferred from the amplitudes noted above (or the discussion in Appendix B) that the 
standard operators by themselves lead to a branching ratio for the fi + K° mode in the range 
B(/j l + K°) std ~ 1 to 10%, while the new operators (related to neutrino masses) by themselves 
can lead to B(fi + K°) new from a few up to 40 or 50%.Q 

With contributions from both the standard and the new operators present, and compa- 
rable in magnitude, one must of course add the two contributions allowing for interference 
between them. Adding the contributions of the standard and the new operators to the am- 
plitudes for all three decay modes (v T K + , T7 fl K + , fx + K°) we estimate the "maximum" and 
the "minimum" of the sum of the two contributions (standard and new) to each of these 
modes, by allowing for the uncertainties in their amplitudes reflected in Eq. (70)-(73). We 
thereby estimate the following range for B(n + K°): 

B(fi + K°) std+ncw « [1 to (50 - 60)%] p (74) 

where p denotes the ratio of the squares of relevant matrix elements for the fi + K° and VK + 
modes. 

If one uses the chiral Lagrangian method as a guide |]58| , |34]| , one would obtain p = 
[\-{m p /m B )(D-F)} 2 x [2/3(m p /m B )Dq+m p / (3m B ) (D+3F)}~ 2 , where D ~ 0.76, F ~ 0.48 
and tub — 1150 MeV, and q = C(usdv) / C{udsv). The entity q denotes the ratio of the 



15 In estimating the upper and the lower ends of the sum {F(V T K + ) + T(l7 fl K + )} and thereby 
the expected range of B(p + K°), we stipulate that the maximum (or minimum) of A(V T K + ) is 
accompanied by a median rather than the maximum (or minimum) value of A(v^K + ) and vice 
versa. This is because the various terms contributing to the two amplitudes are unlikely to be 
entirely constructive (or destructive) simultaneously for both of them. Thus the logical maximum 
and minimum of B(p + K°) lie beyond the stipulated range quoted above. 
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coefficients of the d = 5 amplitudes leading to spinor contractions (us)(du) and (ud)(sis) 
respectively (see Appendix B). In the SO (10) model (developed here), we find typically 
q pa 0.2 to 0.6. Thus the chiral Lagrangain method (CLM) would suggest p ~ 1/4 to 1/5. 
It should however be noted that CLM leads to matrix elements that are inconsistent with 
the results of lattice calculations [4I[] by factors of 1.5 to 4, for proton decaying either into 
(£ + n° and Vir+) or £ + K° or both for all values of (3 H ^ (0.006 GeV 3 )(l/2 to 3) (see the 
discussion following Eq. (116) for the definition of f3jj.) For instance, for (3h — 0.003 GeV 3 , 
the p — > 7T° and p — > tt + amplitudes are larger by a factor pa (1.5 — 2) for CLM relative to the 
lattice calculation, while that for p — > K° is smaller by a similar factor for CLM relative to 



the lattice value Lattice calculations of vK + are apparently not available at present. 
Evidently, one cannot regard the estimate of p based on the chiral Lagrangian method as 
entirely reliable. Even a modest 25 to 30% correction^ to the matrix element of VK + and 
p + K° (compared to CLM values) can alter the estimate of p by more than a factor of 2. 
In the absence - presumably temporary - of a reliable calculation, one should remain open 
to the possibility of larger values, say p ~ 1/2 to 1. Clearly, it is only such larger values of 
p, together with the estimate of B(p + K°) presented in Eq. (74), that permit optimism for 
discovery of the p + K° mode. 

Returning to the estimate of Eq.(74), we find that typically for a large range of parame- 
ters, which correspond to the uncertainties in the amplitudes exhibited in Eq. (72)- (73), the 
branching ratio B(p + K°) can lie in the range of 20 to 30% (if p ~ 1). Thus we see that the 
p + K° mode is likely to be prominent in the SO (10) model presented here, and if p pa 1 it 
can even become a dominant mode. This contrasts sharply with the minimal SU(5) model, 
in which the p + K° is expected to have a branching ratio of only about 10 -3 . In the 5O(10) 
model, the standard operator by itself gives a branching ratio that is at least an order of 
magnitude larger than the SU(5) value, while the potential prominence of the p + K° mode 
arises only through the new operator related to neutrino masses. Thus the p + K° mode of 
the proton decay serves as a signature for the sort of 50(10) model of fermion masses and 
mixings explored here. 

Because of the hierarchical nature of the Yukawa couplings in the Dirac as well as the 
Majorana sectors, proton decay into e + K° and e + ir are highly suppressed in our scenario 
compared to the VK + or the p + K°. 



VII. SUMMARY AND CONCLUDING REMARKS 

One major goal of this paper has been to understand the masses and mixings of the 
neutrinos, suggested by the atmospheric and the solar neutrino anomalies, in conjunction 
with those of the quarks and charged leptons. Adopting familiar ideas of generating eigen- 
values through off-diagonal mixings, we find that the bizarre pattern of masses and mixings 



16 The chiral Lagrangian method has played a role in Ref. [34] and thus in our discussion in 
Appendix and Sec. VI on proton decay rate. We note that corrections of 20 to 30% in the 
amplitude would not however alter the total proton decay rate by more than a factor of 1.5 to 
atmost 2. 
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observed in the charged fermion sector can be adequately described (with ~ 10% accuracy) 
within an economical 50 (10) framework. A concrete proposal was presented that provides 
five successful predictions for the masses and mixings in the quark and the charged lepton 
systems. The same description goes extremely well with a value of m Vr ~ 1/20 eV as well as 
with a large v ll — v T oscillation angle (sin 2 26>p c ~ 0.82 — 0.96), despite highly non-degenerate 
masses for the light neutrinos. Both these features are in good agreement with the SuperK 
results on atmospheric neutrinos. 

The other major goal of this paper has been to revisit the previously noted link between 
neutrino masses and nucleon decay |T2] in the light of the SuperKamiokande result. We 
find that the mass of v T (~ 1/20 eV), together with the large — v T oscillation angle 
implied by the SuperK result, suggest a significant enhancement in the standard as well 
as in the new (neutrino mass related) d = 5 proton decay operators, relative to previous 
estimates, including those of Ref. [|l2j. There are many uncertainties in the prediction for the 
proton decay rate, including ones arising from uncertainties in the SUSY spectrum, from the 
hadronic matrix elements, from the relative phases of the many different contributions (see 
the discussions in Appendix B), and from the allowed extent of unification scale threshold 
corrections to a^{mz)- Nevertheless, following a detailed analysis of threshold corrections for 
the minimal Higgs system (45h, 16h, 16#) used in our work (see Sec. VI D and Appendix 
D) we found that the standard operator contributions severely constrain the underlying 
model. 

Specifically, we found that for MSSM embedded in 5*0(10), the standard operators, with 
generous allowance for the uncertainties, lead to lifetime estimate T _1 (p — > uK + )^ SM < 
[1.5, 7] x 10 33 yrs corresponding to cases (I, II) (see Sec. V, Eq. (37) for the origin of these two 
cases). In the process, by combining constraints from the observed limit on proton lifetime 
together with a reasonable upper limit on threshold corrections, we found that tan (3 must 
be rather small (< 3,8) for cases (I, II), in MSSM. For larger values of tan/? ~ 20, one 
must turn to the Extended Supersymmetric Standard Model (ESSM) embedded in 5*0(10). 
This allows for two extra (vector-like) families at the TeV scale, and has been motivated on 
other grounds [46]. In this framework, the standard d = 5 operators still represent perfectly 
viable sources for proton decay. We find typically T _1 (p — > uK + )^^ M < (1 — 5) x 10 33 yrs 
for tan/3 > 20 (Sec. VI.D). 



As observed in our earlier work fll2 |, one very important consequence of quark-lepton 
unification is the likely existence of new d = 5 proton decay operators that are related to 
neutrino masses. Here we have shown that the "observed" mass v T ~ 1/20 eV and large 

— v T oscillation angle (which go well with theoretical expectations in 5O(10)) enhance 
these operators relative to our previous estimates |T2]. As a result, the standard and the 
new operators appear to make comparable contributions (see remarks following Eq. (64)). 
We have remarked in Sec. VI. F that in some string-inspired models leading to O224 or G2113 
symmetries (rather than intact 5O(10)), the color triplets related to electroweak doublets 
get projected out of the spectrum altogether, and thus the standard d = 5 operators do not 
contribute. In that case, the new d = 5 operators related to neutrino masses will survive 
and dominate. We found (Sec. VI. E) that (given the SuperK result) these new operators 
by themselves lead to proton lifetime T _1 (p — > T7K + ) new with a "central value" of about 
3 x 10 31 yrs, and a range that is perfectly compatible with the observed limit. Assuming 
rather generous uncertainties (i.e., those in (m^/m^), nig, the matrix element /3h), and the 
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amplitude, the new operators by themselves lead to a proton lifetime < (1 — 6) x 10 33 yrs. 
Thus for the MSSM (or ESSM) embedded in 5*0(10), we expect the proton lifetime to be 
shorter than about 10 34 yrs. 

A distinctive feature of the SO (10) framework discussed here is the potential prominence 
of the charged lepton mode p — > fi + K°. In minimal 577(5), this mode becomes prominent 
only for tan f3 > 20, when gluino dressing becomes significant. But such large values of tan j3 
are highly disfavored in SU(5), owing to observed limit on proton lifetime. Furthermore, 
for such large values of tan/3, the Higgsino dressing of the baryon number violating RRRR 
operator becomes dominant, leading anyhow to dominance of the neutrino mode. Thus in 
conventional SU (5)— like supersymmetric unified models, charged lepton decay of the proton 
is (even relatively!) scarce. In the 50(10) framework, on the contrary, the decay p — > fi + K° 
competes favorably with the neutrino mode. This becomes possible primarily because of 
enhancement of the new neutrino-mass related d = 5 operators, in which the amplitudes 
for fi + K° and VK + modes are comparable. Thus observation of the fx + K° decay mode of 
the proton would be very encouraging for the circle of ideas discussed here. Owing to the 
hierarchical Yukawa couplings of the Majorana neutrinos suggested by the solar and the 
atmospheric neutrino data, the modes p — > e + K° and p — > e + 7r° are predicted to be highly 
suppressed relative the p — > VK + . 

While we focused on a specific 5*0(10) example, several of our results are likely to be 
more general. In any scenario with just three light neutrinos, a simultaneous resolution of 
the solar and atmospheric neutrino anomalies argues for hierarchical neutrino masses. In a 
framework that also unifies quarks and leptons (with modest mixing angles), a simple way 
to generate the near-maximal neutrino oscillation angle needed for the atmospheric neutrino 
anomaly is to attribute it partly to the (ji — r) sector and partly to the (v^ — u T ) sector. 
The precise way this division is made is model dependent. As one goes beyond minimal 
supersymmetric SU (5) to a unified framework where small neutrino mass is a compelling 
feature, the predicted lifetime of the proton tends to decrease. In such unified models there 
are various factors contributing to shorten the lifetime (see Appendix C) including (i) an 
enhanced coupling of the muon to the color triplets (relative to the coupling of the strange 
quark), (ii) an enhanced up-quark coupling to the color triplet which scales as ^m c m t rather 
than m c , and (iii) the presence of new operators related to neutrino mass. 

Proton decay has been anticipated for quite some time in the context of unified theories. 
Recent data from SuperK on neutrino mass makes the case for observable proton decay still 
more compelling. With improved searches, especially for the VK + and fx + K° modes, either 
proton decay will be revealed, or some promising and otherwise remarkably successful ideas 
on unification will be called into question seriously. 
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APPENDIX A 



Numerical values of various matrices 



Here we give the numerical values of the matrices V u , Vg and V = V^Vi, where V u 
diagonalises the up quark matrix U of Eq. (28): = V u u^ with (g) and (m) denoting the 
gauge and mass eigenstates, diagonalises the charged lepton mass matrix. We also give the 
numerical values of the matrices (H, A, G, F) as well as their matrix product with V, these 
are relevant for proton decay amplitude calculations. From the fit to the fermion masses 
discussed in Sec. IV- V, we have determined the (approximate) values of the parameters 
(a, T), e, e', 77'; ). For numerical purposes we shall choose their "central values": 



a = -0.1096 r] cb , 7] = -0.1507 rj cb , e = 0.0954 r] £ , 

4.14 x 10" 3 r] rf . 



e' = 1.76 x 10~ 4 r] e ,, 7]' 



With these values, the matrix V u is given by: 



V ~ 



V 



e 

e—cr 



e'(-e + a)\ 
e + a 

~(e + a) 1 / 



1 



(75) 



/ 1 -0.061^/ 

0.061^/ 1 
\8.6x 10-Vb^' 0.0U2 Vcb 



-3.61 x 10~ 5 ^ e A 
-0.0142t7 c6 

1 / 



An analogous expression for Ve is obtained by the replacement a — > i], e — > — 3e and e' 
in the first part of Eq. (76). The numerical values of Ve and V = V^Ve are: 



(76) 
-3e' 



V, 



( 1 0.07?^ 5.6 x l0-% b ri ri , \ 

-0.07^/ 0.905 -0.4369^ 
^ -0.031^77^ 0.43697^ 0.905 / 



1 0.0527?,/ + 0.0777,,/ -0.02677 cb 77 e / 

V' ~ I -(0.052?7 e / + 0.07^) 0.911 -0.4227r7 cb 

-0.030?7 c6 77^ 0.4227r7 cfe 0.911 



Numerical values of the matrices (H, A, G, F) and their products with V are: 
H~h 33 



-1.1 x 10~ 5 
-1.8 x 10"V 
-5.8 x 10~ 3 7? c6 77 e 



-1.8 x 10-% 
-2.9 x 10~ 3 
-0.0954r7 cfe 



-5.8 x 10' 3 r] cb r] e 
-0.0954?7 c6 
1 



(77) 



HV ~ h 



33 



/ 1.8 x 10~Vw 

3.0 x 10"%/ 
\2.9 x lO^r]^ - 0.024t7 c6 77^ 



-2.6 x 10-% 
-0.043 
0.3273r7 cb 



-5.8 x 10" 3 i] cb r] t r \ 
-0.0954?7 cfe 
1 



(78) 
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3.1 x 1(T 5 3.4 x 10~%/ 5.8 x lO -3 ^^ 
A r: /,,,. j 3.4 x 10"%, 2.7 x 10~ 3 0.0954^ cfe 
5.8 x 10~% b ?7 e / 0.0954r? cfe -2.7 x lO" 3 



AV ~ /i 



33 



/ -1.9 x 10'Vw + 3 - 5 x 10_5 
3.0 x 10~%/ + 5.6 x 10- V 
V -6.7 x lQ-% b ?v + 8.6 x KTSfc^e 



2.8 x 10"%/ 5.8 x 10-%^ 
0.043 0.0954?? cfe 
0.0954?7 cfe -0.043 



, (79) 



tan7G ~ h 



33 



5.4 x 10'Vw 41 x 

4.1 x 10-%/ -1.2 x 10- 3 
-2.5 x 10- 3 r] cb r] e , -0.0411^ 



-2.5 x 10 _ % b ?7 e / 
-0.041 lr^ 
1.2 x 10- 3 



t&wyGV' ~ h 33 



/ {3.7 x 10 _ V»V {-1.1 x 10~V {- 2 - 5 x lO^-Vb^ \ 

-2.9 x 10- 4 } +4.1 x 10-%/} -1.8 x 10-%6^/} 

5.4 x 10-%/ -0.019 -0.0411?7 cf) 

V -3.7 x 10"%^/ + 2.9 x 10-% 6 ^/ -0.0411?7 cf) 0.019 J 
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/x + 7.3x 10- 7 + 1.0 x 10~% b y 1.7 x l(TVw - 0.06^/x 0.06^/y \ 

1.7 x 10~%/?7 cb y - 0.06?7 e /x 0.0284^ y + 0.0U2r] cb 

V 0.06^/y y + 0.0142^ 1 / 



(81) 



APPENDIX B 

B.l. Numerical evaluation of the full amplitudes for proton decay 

In converting the superpotential given in Eq. (41) into a proton decay amplitude, the 
first step is to dress two of the sfermions among the four superfields in each term of Eq. (41) 
i.e., to convert them into fermions. The dominant contribution to the amplitude arises from 
the dressing involving the charged wino. Here we list the 12 terms that arise in dressing 
each one of the 7 terms of Eq. (41) with a final state neutrino. The four fermion amplitude 
after wino dressing of a generic superpotential term 

W = M~g (u T Fd') {u T G£ - d'TGu'} (82) 

is 

A(p - VX) = M£e afh x [F n G 2i KO( S '%') (f(c, I) + f(u, d)) 

- F 22 G u {u a ^){s'^{) (f(c,l) + f(c,d)) - F l2 G n {u a d'^v[) (f(c,l)+f(u,d)) 
+ F 12 G 2l (u a s , ^i/ l ) (f(c,l) + f(d,c)) - F 13 G ll (u a d'f)(b'^' l ) (f(t,l) + f(u,d)) 
+ F l3 G 2l {u a b'^u[) (f(c,l)+f(d,t)) + F 13 G 3l (u a b'^(b'^ l ) (f(t,l) + f(d,t)) 

- FnGui^b" 3 )^) (f(c, I) + f(t, d)) - F^G^s'^vl) (f(c, d) + f(t, I)) 

- F^Guiu^b'^vl) {fit, I) + /(*, d)) + F^G^d'^Wul) (f(t, I) + f(u, d)) 

+ F l2 G 3l {u a s'^v[) (f(t,l) + f(d,c))} (g) . (83) 



40 



The dressing functions f(a,b) are defined later in Eq. (93). In writing Eq. (83), use has 
been made of the symmetric nature of F (F^ = Fji). As described in the text, the notation 
d' — VcKud and v' = Vq KM v has been adopted. 

A simpler expression can be obtained for proton decay amplitude with a charged lepton 
in the final state (for the same generic superpotential term as in Eq. (82)): 

A(p -> tX) = M£e a(h x [{u a d p ){uUl)[V cd {G n F 21 - G 2l F n ) + V td (G u F 31 - G 3l F u ) 

+ (u a s^£-)[V cs (GuF 21 - G 2l F u ) + V ts (G u F 31 - G 3l F u )} . (84) 

The amplitude with a hat is the full four-Fermion amplitude divided by the loop function 
(2/), defined in Eq. (110) and (111) of Appendix C. Note that the down type quark fields 
appearing above (without primes) are the physical ones. 

With the numerical values of the Yukawa coupling matrices relevant for the color triplet 
exchange, we can use Eqs. (83) and (84) to compute the decay amplitude for any given 
channel. Since there are seven terms in Eq. (83), for the neutrino mode, there will be a 
total of 7 x 12 = 84 terms to be summed. This is most efficiently done numerically using 
Mathematica. We now present some of the dominant amplitudes. For these estimates we 
define V cd = 0.22r] cd ,V td = 0.006r] td ,V ts = 0.04i] ts and V ud = V cs = 1. The are the 
unknown phase factors, they can also be used to vary the central values of the CKM mixing 
angles adopted. We drop terms which are smaller by more than an order of magnitude 
compared to the leading term in each set. 

A\{ud){sv T )\ ~ M^h 2 33 [1.9 x 10- 5 r] cd r] ts r) £/ - 9.0 x 10~V^ C ^' 

- 6.1 x l0-% d ri ts ri cb ri e , + 4.8 x 10" V - 2.5 x 10~V 

+ 2.1 x 10'%^ + 3.0 x r] cd r] ts r] rf + 2.2 x 10~ e r) cd r) cb r) vl } 

+ (Mi 6 tan7)~ 1 /i 33 /33[3.1 x 10~ 7 r) cd r) ts r) e , + 6.0 x 10~ 7 r)t d r) ts r) cb r) e/ 

+ 4.3 x 10~ 7 7itdrhsVd>Vn' + 1-5 x 10 _ Vi?V 

+ 2.3 x 10 -7 'rjcdVtsVv' ~ 0.041 l^x + 3.1 x lO" 5 ?^^ 

+ 2.2 x I0" 5 r] cd r]t s r] cb r] e >y + 1.1 x lO -5 ^^?/ 

+ 1.6 x 10~ 5 rj cd r) ts r} cb r} v <y + 1.1 x 10~ 5 r] cd r] vl y 

+ z{-2A7 x 10" 4 ^ cb - 1.28 x l(rVi 

+ 9.94 x 10- 5 ms r] cb r] e/ + 7.23 x lO -5 ?^^}] . (85) 

A[{ud){dv T )\ ~ M^h^[2.9 x 10-% dVtdVe , - 2.0 x 10~V^] 

+ (M 16 tan 7 )- 1 /3 3 / i 33[4.7 x KTVi^' + 8.9 x 10~V^ 

+ 1.8 x l^r] cb 7] rf + 1.3 x 10~% d + 6.5 x 10- s r] cb r] vf 

+ 6.8 x lO'Vd^W - 1-8 x 10~ 8 r/ c6 ?7 r? / - 9.0 x 10~ 3 r] cd r] cb x 

- 6.8 x 10~ + 3.3 x l0- 6 rj cd rj td rj cb rj e/ y 

+ 4.8 x 10^r] cd r] td r] cb r] v ,y + 2.5 x 10~\>y 

+ z{-hA2 x lOScdritdVcb - 2.82 x 10~ 5 + 1.49 x 10~ 5 ^ c& ^ 

+ 1.08 X VtdVcbVv'}} (86) 
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A[{ud)(su^)} ~ x 10-% d r) ts r) £ , + 8.5 x lO^V^^e' 

- 3.4 x 10"V^' + 2.1 x l(TV»7c& - 4.7 x 1(T 6 7^WV 

- 7.0 x r] cd r] ts r] cb r] v , + 5.1 x 10' 6 r] cd r] V - 2.1 x 10" Vw 

- 1.9 x 10-% s r] cb r] € ,r] ri ,] + (M 16 tan-f)' 1 f 33 h 33 [2.Q x 10~ 7 rj td rj ts rj e/ 
+ 1.4 x i] cd r] ts r] cb r] e > - 9.9 x lO -7 ^^/ 

- 3.5 x 10~ 7 rjurjarj^ - 5.2 x 10~ 7 rj cd rj ts rj cb rj v/ - 1.8 x 10~ 7 r] cd r] v , 

- 1.4 x 10~ 7 rjtsrjchrje'rjr,, - 1.6 x 10~ 3 r) ts r] cb x - 1.8 x 10~ 2 a; 

+ 9.2 x 10- 6 r) cd r) ts r) e ,y - 1.4 x 10^r] cd r] cb 7] t ,y + 2.5 x 10~Vy 

- 2.5 x 10~ 5 r] td ri v iy - 3.6 x 10~ 5 rj^rjisTj^y + 2.6 x 10~ 5 r) cd r) cb r) v 'y 

- 1.0 x 10~ 5 11^^^ - 7.1 x 10~Vf>W^'2/ 
+ ^{-9.86 x KTV^** - 1.11 x 10~V 

- 5.79 x 10- 5 Mc6 + 4.20 x lO" 5 ?^}] (87) 

Similarly, the amplitudes A[(«s)(rfi/ T )], A[(ws)(<ii/ M )] and A[(-u<f)(<iz/ M )] can be computed. 
We do not display these results here, since the amplitudes for these operators are always 
somewhat smaller than the ones displayed. Furthermore, the matrix element for the former 
two turn out to be suppressed by about a factor of 3. 

Turning now to the charged lepton decay amplitude, 

A[(us)(ufi-)]) ~ AQ 1 ^[V r „{-8.4 x 10~ 7 - 2.1 x 10"V%'} 
+ V^{5.2 x 10" 5 -4.8 x 10-\ bVe , W } 
+ (M 16 tan 7 )- 1 / 3 3^33[{K ;s {-5.0 x 10" Vw + 1.8 x 10~ 2 x 

- 7.1 x 10- 6 yr] cb ri v ,ri e ,} + V ts {9.3 x lO" 7 ^ - 3.6 x 10-% b r] v 'Ve' 

- 0.0411^ + 6.4 x 10~ 5 y + 2.5 x lO -4 ^^/} 
+ z{4.2 x 10 _5 ?7^77 e , - 1.59 x 10~%r] e > 

- 1.66 x 10-%^, - 5.88 x lO -5 ^}] . (88) 

A[(us)(ue-)\ ~ hl z M^;\V cs {SA x 10~%, + 1.5 x 10"V} 

+ \4 S {8.3 x 10- 7 We , - 6.9 x 10- 6 77 c6 77^ + 3.3 x 10~V^} 

+ (M 16 tan 7 )- 1 / 3 3^33[K ;s {3.5 x 10" V - 1-7 x 10" 

+ 5.3 x 10~ 3 ?v:r + 5.0 x 10"V 6 ^y} + V te {8.3 x 10~\ b r] e/ 

+ 3.1 x 10^ 7 r] cb i] v r + 2.5 x 10~ 7 r) cb r) e ' + 3.7 x 10' 4 r] cb i] v rx 

+ 2.9 x KTVtftyz + 2.2 x 10" 5 ^/y + 1.8 x 10 _ Vy} 

+ Z {-1M x lO'V^'W}] • ( 89 ) 

The decay amplitude for p — > 7r°/i + can be obtained from Eq. (88) by the replacement 
V cs — > and Vt s — > 14^- The amplitude for p — > e + 7r° may be obtained from Eq. (89) by 
the replacement: V cs -> V^, -> Kd- 
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B.2 Representative contributions to proton decay amplitudes and their estimated 

magnitudes: 



In this subsection, we exhibit in detail the evaluation of a few representative (dominant) 
contributions to proton decay amplitudes. The full set of contributions are listed in the 
preceding subsection. We also estimate the magnitudes of the full amplitudes allowing for 
uncertainty in the relative phase of the different contributions. 

To obtain some of the leading terms in the results exhibited above, consider the part 
containing the v' field - in the first term of (Eq. (41)), which we specify by subscript 
Iv'. 



W$ = -M eS e aP7 H ij (HV% l (u?df)(d%u' l ) 



(90) 



Here (a, (3, 7) denote color indices. The transpose symbols on the relevant fields are dropped 
henceforth. We will first consider case I discussed in Sec. V for which e' 7^ 0. The first 
and the second terms in Eq. (82) arise by making the choices (A and B) and (C and D) 
respectively, for the indices as given below: 



(ij)(kl) 



(21) (23) 

(22) (13) 

(23) (13) 
(21)(33) 



(Choice A) 
(Choice B) 
(Choice C) 
(Choice D) 
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Note that A and B are related by the interchange of j <-> k, and similarly C and D. 
Choice A: (ij)(kl) = (21)(23): Starting with the operator (u^df)^'^^) in the superpo- 
tential, two fields need to be utilized for wino dressing. Which ones remain external gets 
determined as follows. The field d'2 must remain external, yielding a strange quark s 7 , 
accompanied by the CKM factor V cs ~ lf£] (If d' 2 were dressed instead, it would yield an 
external charm quark, which is kinematically disallowed.) Thus both an d df must be 
utilized for dressing, which, after conversion of c and d[ at the wino vertex, respectively, 
yield V c dd a = ij] cd c )d a and V^u 13 , where 6 C — 0.22. As discussed before, r\ cd = ±1. r] ud 
and f] us are chosen to be +1 by convention. Thus, the corresponding contribution to the 
four-fermion proton decay operator leading to V 3 emission is given by 



A 0) 



-f(c, d) 



h 2 ' 
M 



cff 



(OcVcd)ea^H 21 (HV') 23 (d a u%s~<is 3 ) 
/(c,d)(3.8 x lQ- 6 )r lcd rie>ricbe a p~ / (d a u> 3 )(s"<v 3 ) 



(92) 



where 



17 If d' 2 is used to yield an external d~quark, and U2 (after dressing) an s-quark, one would obtain 
the four-fermion operator (su){dv 3 ). It turns out that the matrix element of this operator is 
suppressed by about a factor of 2-3 compared to that of {du){su 3 ). See [|||]. 
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\4n J ml - m 2 b \ml - m 2 w m 2 w J Kin J 

In getting Eq. (92), we have used the numerical values of the elements of H and HV given 
in Appendix A (Eq. (78)). The phase factors rj e i and r] cb are ±1. Likewise, confining still to 
the index combination A, the contributions from the remaining five standard operators of 
labeled by the subscripts (II- VI) are found to be: 

(Afly, Afl iy , Aty^AW, A$„)~(-3, +2, -6, -0.4, -0.8)4? (94) 

Combining the contributions (92) and (94), the total contribution of the terms containing 
v' in the first six operators (I to VI) in W L is given by 

77^^(2.74 x 10- 5 )e^ T (d a ^)( S V 3 ) . (95) 

Next consider choice B. In this case, since one starts with the operator eap-yiu^d'^idi V3) , 
using the dressing as in case A, one gets an extra minus sign compared to A, owing to the 
color factor. Following the same procedure as above, one obtains 

VcdVe'Vcb(-l.S2 x 10' 5 )e a ^(d a u^u 3 ). (96) 

Adding Eq. (95) and (96), the total contribution of terms containing v' in the first six 
operators of for the index combinations A and B, leading to V 3 K + emission, is given 
by: 

(A (1) + B (1) ) /to viy- 

This checks with the second term in Eq. (85) obtained numerically, up to an overall sign, 
which is due to the difference in the sign of Eq. (96) compared to Eq. (41). Likewise, the 
corresponding contributions for the index combinations C and D from the first six operators 
in is found to be: 

(C-W +£>«), to VI ~ 

This is the first term of Eq. (85) - up to an overall sign. Note that Eqs. (97)- (98) add 
constructively because m s = —r]^. We still have to include (i) the contribution that arises 
by interchanging % <-> j while keeping (k, I) fixed (see Eq. (83)) and (ii) that from the terms 
containing charged lepton field I in W^ L \ The latter contribute to neutrino emission by 
dressing the charged slepton fields (£). Including contributions from both (i) and (ii), one 
can verify that the net contribution of each kind - such as A, B, C, D listed above - is 
obtained by simply making the following substitution in Eq. (97)- (98), etc: 

/(c, d) - /(c, d) + f(c, I) . (99) 



(i) 



to viy 



h hf(c, d) 



M, 



cS 



D i to viy 



hhf{c, d) 



M, 



eff 



h U (c, d) 



eff 



VcdVe'Vcb(0.9 x lO-^e^cfV 3 )^) . (97) 



hhf(c, d) 



M, 



eff' 



VcdVe'VU-l-9 x 10^ 5 )e a/37 (^)( S V 3 ) 



(98) 
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In arriving at this result, use has been made of the symmetric nature of the matrices 
H, A,G, F, as well as color antisymmetry. Thus, adding Eq. (97)-(98), and making the 
above substitutions, the net contribution to V 3 K + decay mode from the first six terms of 
W^ L \ for the index combination [A+B+C+D + (i «-> j)] is given by: 



[A + B + C + D + (i++j)] I to VI 



h 2 

' l 33 



f(c,d) + f{c,£) 



VcdVe'Vcb X 



(2.8 x 10- & )e a/37 (^)(s 7 z/ 3 ) . (100) 



Having shown how the first two terms of the first bracket of Eq. (85) arise, let us now 
consider all the other terms in this bracket, which are proportional to and estimate 
their net magnitude. It is clear that the joint contribution of the first two terms (Eq. (85)) 
dominates over the other terms in the same bracket by factors of 5 to 10. Although the 
relative signs (phases) of these subleading terms are uncertain (since we do not know the rjij), 
we can see that their combined magnitude even if they all add constructively, is less than 
about (1.5 x 10 -5 ), which is about half of that given by the first two terms of Eq. (85). Not 
knowing the relative signs (phases) we conclude that the net contribution from the standard 
operator (first six terms of W) is given by:0 h 2 33 [2f(c,d)/M eS }(2.8 x HT 5 )(l/2 to 3/2). 

Turning to case II (e' = 0,5 ^ 0, see Sec. V), the net contribution from the standard 
operator to the same amplitude is obtained simply by setting T] e i = in the first square 
bracket in Eq. (85). Since 5 ~ 10~ 5 ~ e'/20, contributions proportional to 5 can safely be 
neglected compared to the other terms. Noting that r] ts = —rj cb , the said contribution is thus 
given by [(4.8 + 2.5) x 10~V + (3 + 2.2) x l(rV»7«i»ty] = (0.73 + 0.55?^)^ x 10~ 5 . 
Allowing for unknown relative phase as above, this contribution will vary in magnitude 
between (1.28 to 0.18))xl0~ 5 , which may be represented approximately by (0.7 x 10~ 5 )(1 ± 
0.7). 

Thus we estimate that the contribution of the standard operator to the V T K + amplitude 
for the cases I and II are given by: 



A(u T K + ) std 



hl 3 2f(c, d) 



M, 



cir 



,„ b\, 7 \ f 2.8 x 10~ 5 (± to |) -Case I 1 f . 

^(d U^U Z ) X( 07x 1Q _ 5( ^ ± Q _2j _ ^ n | . (101) 



Note that the standard amplitude for case II (e' = 0) is typically smaller than that for case 
I (e' 7^ 0) by a factor of 3.5 to 4. Of course, if there is a near cancellation for case II for the 
V T K + mode, so that the curly bracket in Eq. (101) (for case II) is as small as 0.2 x 10~ 5 , 
then contributions from the new operator (related to neutrino masses) to V T K + as well as 
contributions of the standard and the new operators to the V IJL K + mode would turn out to 
be far more important, as discussed below. 

Contributions to the V^K + mode from the standard operator (First six terms in 
W) 



18 The reader will note that these unknown phases (or signs) are just the analogs of those which 
enter into t versus c contributions in minimal SU(5), which is represented by a parameter ytx etc, 



see Ref. [HJ3§ 
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Contributions to the V ll K + mode from the standard operator for the cases I and II can 
be estimated by using Eq. (87) in a manner similar to that noted above. Collecting terms 
proportional to i] e r and those independent of it from the first square bracket of Eq. (87), 
and using r) ts = —i] cb = ±1, the said contribution is given by: [{—1.2r) cd — 0.35r] td r] ts — 
(0.21 - 0.19)^/}^/ + {1.2r] cd r] v i - 0.21 - 0A7i] td r] cb r] v r}] x 10~ 5 . Allowing for uncertainty of 
relative phases, a reasonable estimate of the magnitude of this contribution is thus given by 
[1.5(1 ± 0.5) x 10~ 5 ] for case I (e' ^ 0) and by [1.2(1 ± 0.5) x lO" 5 ] for case II (e' = 0). We 
thus obtain: 



^4(V^ + )U- 



hh2f(c,d) 



M, 



cff 



\\ ) (1 ± 0.5) x 10-^(^X5%) . (102) 



The upper and lower entries refer to cases I and II respectively. We observe that the 
amplitudes for V^K + and V T K + modes have very similar contributions from the standard 
operators; the former is expected to be leading for case II and the latter for case I. 

Combining the rates for the two modes, we would expect that the total proton decay 
rate, represented to a good approximation by the sum \T(V T K + + T(v lJL K + )] ) would typically 
be proportional (for central values of the amplitudes) to [(2.8) 2 + (1.5) 2 ] ~ 10 for case I, 
and to [(0.7) 2 + (1.2) 2 ] pa 2 for case II. In short, in so far as contributions from the standard 
operator is concerned, proton decay rate would typically be suppressed by about a factor of 
5 for case II compared to case I. 

T((v T K+) atd + T((v^K + ) std ] Casc ii _ 1 
V{{v T K+) std + V(iv ll K+) std } c&sc j ~ 5 ' 

This is an interesting and significant effect on proton decay, having its origin entirely in the 
fermion mass matrix. 

Contribution to V T K + and V^K + mode from the new operator related to neutrino 
masses (the seventh term in W^): 

Confining first to the part containing v' in the seventh term of (see Eq. (41)), we 
have: 

W$is = -^UGV%{utdf{d'M) . (103) 

Since F 33 is the leading element of F, first consider the following combination of indices: 
Choice E: (ij)(kl) = (33) (13): Starting with the operator (u 3 df)(d'iU 3 ), it is clear that df 
must be external. This yields Vts-s 13 (if d 3 is dressed, it would lead to external top, which 
is forbidden); thus and dj must be dressed to respectively yield V td d a and u 1 . Thus one 
finds: 

E W „ F ^p^ Vtd V ts f{t,d)e aM {d^){s^,) 

M 16 



/33 ^33 



-(QVcbVe' + ^VcbVv^VtdVts x 10 7 



X 



\Miq tan7^ 

CKMf(t,d)e aPl (d a u^iy 3 ) . (104) 

Here, we have used values of F^ and (GV')ki tan 7 listed in Appendix A (Eq. (80)- 
(81)), and have put V td = 0.00Qr] td , V ts = 0.0Arj ts . The factor ( KM = 1/2 to 1 denote the 
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combined uncertainty in the two CKM elements. Note that the two terms in Eq. (104) 
check with the second and the fourth terms in Eq. (85) proportional to (M 16 tan7)~ 1 . 
Some of the other leading contributions arise by choosing (ij) = (23) or (32), varying (kl). 
These are proportional to y and is given by the seventh through the eleventh terms in the 
second bracket of Eq. (85). Note that the tenth and the eleventh terms subtract (because 
VtsVcb — — 1) ; and thus have a combined magnitude ~ 0.5 x 10~ 5 y, while the seventh and the 
eighth terms subtract to have a net magnitude of 0.9 x 10~ 5 y. Thus, allowing for unknown 
phases of each of these contributions and also of the 10th term, we estimate that the total 
contribution from the y dependent terms is (1.3 ± 1) x 10~ 5 y ~ (0.7 ± 0.5) x 10~ 6 , where 
we have put y ~ 1/20 (see Eq. (25)). Now the y-independent terms in the second square 
bracket of Eq. (85), including that exhibited in Eq. (104)), are each individually smaller 
by factors of 1.5 to 2.5 compared to 0.7 x 10~ 6 except that the term proportional to x is 
nearly 4 x 10~ 6 (l/5 to 1) for x ~ 10 _4 (l/5 to 1). Although the relative phases (signs) of the 
various contributions, which depend on rjij, are not known, allowing for possible cancellation 
between some of the smaller terms, it seems most plausible (as a conservative estimate) that 
the second square bracket in Eq. (85) has a net magnitude m (1.5 x 10~ 6 )(l/2.5 to 2.5) for 
case I. Estimating similarly, one obtains a very similar magnitude for case II as well. Now, 
including the contributions from interchange of i <-> j, keeping (k,l) fixed, as appropriate, 
and that involving charged lepton dressing, as before, the net contribution from the new 
operator (the seventh term in W^), related to the neutrino masses is thus given by: 

A(u T K + ) QCW ~ f j 33 f 33 ) [(1.5 x HT 6 )(l/2.5 to 2.5)1 [f(t, d) + f(t, I)] 
\Mi 6 tan7y L J 

e a/3 ,(d a u^u 3 ) (105) 

for case I as well as case II. 

Using Eq. (87), the contribution of the new operator to the V^K + mode for both cases 
I and II is estimated to be: 

A{V,K + ) acw « ( J 33h3 ) [(3 x 10~ 6 )(l/2 to 2) 
\A2i 6 tan7/ L 

x[f(t,d) +f(t,l)}e a f 3j (d a u^)(s^ fl ) (106) 

Here we have used x ~ 10~ 4 (l/5 — 1) and z ~ (1/400 — 1/200), in accord with discussions 
in Sec. V. Thus, for the new operator, the V /Ji K + mode would typically supercede the V T K + 
mode, and the amplitudes for either mode remain essentially the same for both cases I and II. 
It is worth noting that the main source of the relative enhancement of the v fl K + amplitude 
is due to the fact that the dominant y-independent terms add for V I1 K + (see Eq. 88), and 
they subtract for V T K + (see Eq. (85)), if one uses r] ts = —r] cb . 
Charged lepton decay modes: 

The amplitude for proton decay into charged leptons, p — > £ + X, is given in Eq. (84) for 
a generic superpotential of the form Eq. (82). The number of terms are fewer compared to 
the amplitude for neutrino mode, so they can be summed easily. The full result is given in 
Eqs. (88)-(89) for the muonic and electronic modes. We will briefly describe the origin of 
the dominant terms in the amplitude involving charged muons (Eq. (88)). 
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The standard d = 5 amplitude receives two contributions, one from the exchange of c — d' 
squarks (the term proportional to V cs in Eq. (88)), and one from the exchange of td' (the 
V ts term). They are comparable numerically. Note that unlike in minimal SU(5) models, 
where these contributions undergo a unitarity cancellation, here the amplitude is not much 
suppressed compared to the neutrino mode. 

Following the same procedure as discussed in the other cases, the contribution of the 
standard operator to the n + K° mode (given in Eq. (88)) is estimated to be 



A(fi + K°) std ~ 



ht 3 2f(a,b) 



M, 



off 



3 x 10- 6 (l/2 to 2) - Case I 
(2 ± 0.8) x lO" 6 - Case II 



(107) 



Thus the standard amplitudes for the [i + K° mode are nearly the same for both cases I and 
II. They are typically smaller however than those for both V T K + and V fl K + modes by a 
factor of 3 to 10 (compare Eq. (107) with Eq. (101) and (102)): 



A^ + K\ td /A{V,K- 



Jstd 



1/3-1/10, e = n,r . 



(108) 



We see that for the standard operators, the charged lepton (/i + ) mode, though relatively 
enhanced compared to minimal SU(5) is still expected to be suppressed in the rate, even 
in the SO(10) model presented here, by one to two orders of magnitude, compared to the 
neutrino modes. 

The contributions of the new operator (related to neutrino masses) to fi + K° may also 
be estimated using Eq. (88). Suppressing relative phases, it is found to be: 



A{fi + K°) 



33 



M 16 tan7 / 

^33 ^-33 

Mi 6 tan7 



1.8 x 10~ 6 (:r/10- 4 ) + (2 to 10) x 10~ 7 
10" 6 )(l/3 to 2) (for I andll) 



(109) 



In above, we have used z ~ 1/300 (in Step I) and x ~ (1/3 — 1) x 10~ 4 (in Step 2) in accord 
with discussions in Sec. V. Thus, confining to the contributions of only the new operator, 
we would typically expect: A(fi + K°) ncw / A{p^K + ) ncw m 1 — 1/3 for i = fi, r. (Compare Eq. 
(109) with Eq. (105) and (106).) This means that if somehow the standard operators were 
absent (or highly suppressed) and the new d = 5 operators were the only (dominant) source 
of proton decay (a possibility discussed in Sec. VI. F), the charged lepton fi + K° could be 
comparable to the neutrino mode, having a branching ratio of nearly few to even 50%. As 
discussed in the text (Sec. VI), this is also found to be the case by including contributions 
from the standard as well as the new operators. 

The estimates of the proton decay amplitudes arising from the standard and the new 
(neutrino-related) d = 5 operators for the three relevant modes (V T K + ,V^K + , fi + K°) are 
summarized in Table 1. 



Mode 


(Case I) std (e' ^ 0) 


(Case 11)^ (e' = 0) 


(Case I) new (e' ^ 0) 


(Case II) new (e' = 0) 


V T K + 


2.8 x 10" 5 (1 ±0.5) 
1.5 x 10~ 5 (1 ±0.5) 
3 x 10- 6 (l/2 - 2) 


0.7 x 10" 5 (1±0.7) 
1.2 x 10~ 5 (1±0.5) 
(2 ±0.8) x 10~ 6 


1.5 x 10- 6 (l/2.5 - 2.5) 
3 x 10~ 6 (l/2 - 2) 
10" 6 (l/3-2) 


2 x 10~ 6 (l/2.5 - 2.5) 
3 x 10~ 6 (l/2 - 2) 
10- 6 (l/3-2) 
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Table Caption: The magnitudes for the standard and the new operators can be obtained 
by multiplying the respective entries shown in the Table by the factor X = ft| 3 2/(c, d)/M e g 
and Y = f 33 h 33 2f (c, d)/(M 16 tan 7) respectively (see the text for definitions, in particular 
Sec. VI E, where one derives P = f 33 h 33 /(M w tan 7) « 5 x 10- 19 GeV _1 (l/2 to 2), and Eqs. 
(101) and (105) as examples). 



APPENDIX C. CONTRIBUTION TO PROTON DECAY RATE FROM THE 

D = 5 OPERATOR 



In this Appendix we give a general discussion of the numerical value for the d = 5 proton 
decay amplitude. 
Define 

A(j7 i )=A 5 (p^vK + )/(2f) . (110) 

A 5 denotes the strength of the respective four-fermion proton decay amplitudes arising from 
the d = 5 operators in W. One representative contribution to A 5 is given by Eq. (97). The 
net value of A 5 includes contributions from all the d = 5 operators, allowing for all possible 
combinations of the indices {ij){kl). It also includes the wino dressing factor and the relevant 
CKM mixing factors at all vertices. The net values are listed in Appendix B. The quantity 
/ is the average of the two relevant dressing functions. For instance, for the amplitude in 
Eq. (100): 

J = [f(c,d)+f(c,l)]/2. (Ill) 

Assuming approximate degeneracy of the sfermions (for simplicity, see however, remarks 
later), / is the same for all contributions to the amplitude. Then A defined as above is 
just the net strength of the corresponding d = 5 operator in the superpotential, multiplied 
by the relevant CKM factors which occur at the vertices involving the color triplet and the 
wino exchanges. Since / has dimension of (mass) -1 , A has dimension 5. 

For comparison purposes, it is useful to note that for minimal SUSY SU(5) involv- 
ing exchange of color triplet Hq between (cs) and (dv^) pairs, the d — 5 operator has 
strength = (h% 2 hj 2 /M Hc ) ~ (m c m s sin c /(v u v d ))/M Hc , ~ (m c m s sin 9 c /vl) (tan (3/M Hc ) « 
9 x 10~ 8 (tan/3/M# c ), where tan/3 = v u /va and we have put v u = 174 GeV and the ex- 
trapolated values of the fermion masses at the unification scale - i.e., m c ~ 300 MeV, and 
m s ~ 40 MeV. Multiplying further by an additional factor of sin 0c due to conversion of c 
to d at the wino vertex, for the case of SU (5), 

4j(S77(5)) ~ (1.9 x 10- 8 )(tan/3/M^ c ) , (112) 

corresponding to exchange of the pair (c, d). There is a similar contribution involving the 
exchange of the pair (t, d) with the substitution m c sm6c — > mtVtd- 

In calculating the proton decay rate, we will assume the following spectrum of supersym- 
metric particles as a guide: Squarks are nearly degenerate, with masses m 5 - ^ 1 TeV (1 to 
1.5) and wino is lighter than the squarks - i.e., m^/rriq 1/6 (1/2 to 2). Consistent with 
light gaugino masses, we will furthermore assume that TO1/2 (the common gaugino mass) 
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is small compared to mo (the common scalar mass). In this case, starting with universal 
masses for the scalars at the unification scale, we would expect the slepton masses to be 
nearly degenerate with the squark masses at the electroweak scale.[^] The two /-functions 
that enter into the amplitude are then nearly equal (eg: in Eq. (99)), f(c,d) ~ f(c, I)) and 
/ ~ (m^/m^) (a!2/47r). A SUSY spectrum as described, i.e., <C mq ~ m| ~ 1 TeV is 
needed anyway, a posteriori, because without it SUSY grand unified models based on SU (5) 
or 50(10) are likely to run into conflict with the experimental limits on proton lifetime. A 
spectrum of this type is plausible in several scenarios of SUSY breaking.0 

With the strength of (du)(su T ) operator being given by A 5 = A(2f), as in Eq. (109), the 
inverse decay rate for p —>■ V T K + is given by 



V~ 1 {p -> V T K + ) « (2.2 x 10 



31 > 



.07 



0.006GeV 3 



ft. 



H 



(1/6) 



yrs x 
2 



[m 



w 



/m q 



rriq 



1 TeV 



2 x 10" 24 GeV 
W) 



(113) 

o\(d a y L )ui\ P ,k). 

given that one 



Here (3h denotes the hadronic matrix element defined by /?#«£(&) = e aj g 7 
While the range (3h = (0.003 — 0.03) GeV 3 has been used in the past 
lattice calculations yield (3h = (5.6 ±0.5) x 10~ 3 GeV 3 [(41], we will take as a plausible range: 
Ph — (0.006 GeV 3 ) (1/2 to 2). As stands for the short distance renormalization factor of 
the d = 5 operator. In minimal SU(5) it has a central value of 0.67. Although this factor is 
slightly different in the 50(10) model, we shall adopt A s = 0.67 for the 5O(10) model as 
well. 

Note that the familiar factors that appear in the expression for proton lifetime - i.e., 
M Hc , (1 + y tK ) representing the interference between the t and c contributions and tan/3 
- are all effectively contained in A (77). In fact, the analog of Mh c for SU(5) is given in 
our case by two mass scales: M eff and M 16 , representing contributions of the standard and 
the neutrino mass-related operators, respectively. The analog of (1 + ytx) is reflected by 
the uncertainty in the net value of the many terms in Eq. (85), which depends on their 
unknown relative phases (see discussion in Appendix B, where a reasonable estimate of the 
uncertainty owing to the phases is given (Eq. (101) and (105)). 

In minimal SUSY SU(5), A is proportional to (sin2/?) _1 = 1/2 (tan/? + 1/ tan/3), and 
thus approximately to tan (3 for tan (3 > 3 or so. Corresponding to a realistic treatment of 
fermion masses as in Sec. IV and V, this approximate proportionality to tan/? does not 
however hold for 5O(10). The reason is this: If the fermions acquire masses only through 
the 10h in 5O(10), as is well known, the up and the down quark Yukawa couplings will be 



19 Allowing for sleptons to be lighter than the squarks by a factor of 2 to 3 while keeping the 
squark masses to be about 1 TeV would amount to increasing / by about 50%, compared to its 
value for the case of nearly degenerate masses (rriq ~ nip), and thereby enhancing the rate by a 
factor of 2. This possibility is of course not presently excluded. 

20 For example, models of SUSY breaking based on contributions from a family universal anoma- 
lous U(l) D-term, superposed with subdominant dilaton F-term. contributions |)7]] would lead to 
such a spectrum. 
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equal. This would give the familiar t — b — r — v® unification. By itself, it would also lead to 
a large value of tan f3 = m t /mb ~ 60, and correspondingly to a large enhancement in proton 
decay amplitude. Furthermore, it would also lead to the bad relations: m c /m s = m t /m b 
and Vckm = 1- However, in the presence of additional Higgs multiplets contributing to 
fermion masses, such as 16#, which (a) distinguish between the up and the down sectors 
and (b) correct the bad relation mentioned above (see Sec. IV and V), tan/3 can get lowered 
considerably - for instance to values like 10 - 20. Now, with the contributions from (10#) 
and (16#) d to fermion masses, as in Sec. IV and V, the simple proportionality of A to tan/3 
disappears. In this case it is more useful to write A§ or A simply in terms of the relevant 
Yukawa couplings. That is what we have done in writing W (see Eq. (40)) in terms of the 
products of Yukawa coupling matrices like (H) (HV) etc, and likewise for the amplitude. 

It is instructive to compare a typical amplitude for 5*0(10) to the corresponding one for 
SU(5). This ratio is given by: 

^(^ + )ftd (10) „ hi, 2 x 10-s 

A(v^K+) s st u d i5) ~ M ^ L9 x 10- 8 (tan/3/M Hc ) 

« (m Hc /M cfr )(88)(3/tan/3) . (114) 

We have put ftj 3 ~ 1/4 in going from the first line to the second in Eq. (114). Thus we 
see that M eff has to be (88 — 53) times larger (for tan/3 = 3 to 5) so that the 5*0(10) 
amplitude may be comparable to that of 577(5). In other words, if M eff were equal to 
M Hc , one would have a net enhancement by about two orders of magnitude of the d = 5 
amplitude for proton decay in 5O(10) compared to that of 577(5). This large enhancement 
of the amplitude in 5O(10) - apart from the factor of Mij c /M e g - has come about due to a 
combination of several factors: (i) the large off-diagonal coupling of the up-type quarks with 
the color triplets which scale as ^m c m t in the 5O(10) model in contrast to the diagonal m c 
in minimal SU(5), and (ii) the larger muon Yukawa coupling to color triplets (relative to 
the strange quark Yukawa in 577(5)). 

5O(10), however, has a possible source of suppression of the d — 5 amplitude, because 
of the nature of the doublet-triplet splitting mechanism in it. The suppression would arise 
if the mass M 10 of 10' H is considerably smaller than the scale of A (45 h) = Aa, as discussed 
in Section VIA. In this case, M cff = (Aa) 2 /M 10 can far exceed Aa. Since (45 H ) breaks 
5O(10) to SU(3) C x SU(2) L x I3R x (B — L), one of course naturally expects Aa to be nearly 
equal to or somewhat larger than the unification scale Mjj «2x 10 16 GeV. But M 10 can in 
general be one or even two orders of magnitude smaller than Mjj. For instance, 10^ may 
carry a charge that would forbid a mass term (lO^) 2 at the renormalizable level. Such a 
mass term could still effectively arise through nonrenormalizable operators by utilizing the 
VEVs of certain fields fa which do not conserve the respective charge. In this case, the 
mass term M\q may be suppressed by relevant powers of [(<f>i) /M\ where one may expect 
(fa) /M ~ 1/10. If M10 ~ (1/10 - 1/100) Aa, M cff can exceed M v and thus M Hc by one 
or two orders of magnitude. The significance of such a large value of M cS for coupling 
unification is discussed in the text. 

More precise constraints on M ef j and M 16 tan 7 are obtained directly by using limits on 
proton lifetime. They are discussed in the text in Sec. VI. 
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APPENDIX D: THRESHOLD CORRECTIONS TO a 3 (M z ) DUE TO LOW 
DIMENSIONAL MULTIPLETS OF 50(10) 

Given a\ and a 2 at m z , and the MSSM spectrum, the hypothesis of grand unification 
allows us to predict a 3 {m z ). This prediction, however, receives so-called threshold correc- 
tions owing to mass-splittings, at the unification scale My, between submultiplets belonging 
to complete 50(10) multiplets. Such mass splittings are of course induced when 50(10) 
breaks spontaneously to lower symmetries using Higgs multiplets such as 45h and 16h- 

Denoting the one-loop threshold corrections to a^ 1 (m z ) by — Aj, so that a~ l (m z ) = 
a u~ -(b i /27r)\n(m z /M u ) - A h we obtain: Aj = E Q (-^)/( 27r ) ln(Mu/m a ) where h L = 
(33/5,1,-3) yield the familiar one-loop f3 functions for the evolution of the three gauge 
couplings (i = 1, 2, 3) for the MSSM spectrum, and bf is the contribution to the evolution 
from the arth sub-multiplet with mass m a . Using straightforward algebra, the threshold 
correction to a 3 {m z ) is thus found to be: 

Aa 3 (m z ) = [a 3 (m z )] 2 ^A x - yA 2 + A 3 ) (115) 

For illustration only, consider the contributions to Aa 3 (m z ) due to a splitting between the 
doublet and the triplet members in a 5 + 5 of 51/(5) (or equivalently a 10 of 5O(10)) with 
masses m 2 and m 3 respectively, which are both of order My. In this case, one obtains: 

Ax = (l/27r)(3/5)ln(77i2/M^) + (1/2tt) (2/5)ln(m 3 /M^) 

A 2 = (l/27r)(l)ln(m 2 /M [/ ), A 3 = (l/27r)(l)ln(m 3 /M [/ ) (116) 

Using Eq. (116), one then obtains: 

{Aa 3 (m z )} DT = [a 3 (m z )] 2 (9/7)(l/27r)ln(m 3 /m 2 ) (117) 

Using the doublet-triplet splitting mechanism exhibited in Eq. (39) and the fact that 
the MSSM evolution includes the two light Higgs doublets, one obtains the corresponding 
formula for Aa 3 {m z )DT exhibited in the text (Eq. (53)), where (m 3 /m 2 ) is replaced by 
(M cS cos ^ / Mfj) as explained there. 

The threshold corrections owing to mass-splittings within the gauge multiplet and the 
Higgs multiplet 45h can be computed in the same manner. First consider the gauge multi- 
plet. 

The masses of the gauge particles acquiring superheavy masses are given by: 
M 2 (3, 1,2/3) = V(c 2 + a 2 ), M 2 (3,2,l/6) = g 2 (Ac 2 + a 2 ), M 2 (l,l,±l) = Ag 2 c 2 , 
M 2 (3,2,-5/6) = # 2 a 2 , where the numbers in parentheses denote SU(3) x SU(2) x U(l) 
quantum numbers, and g is the unified gauge coupling. Such a splitting leads to a correc- 
tion in a 3 , which at the electroweak scale is given by: 

Aa 3 (m z ) gaugc = - a3 ^ z)2 [751n(4 + p 2 ) - 1051n(l + p 2 ) + 301np 2 ] (118) 

where p = 2c/ a. Eq. (118) includes contributions from the gauge particles as well as from the 
corresponding gauginos and respective Higgsinos which together make 4-component Dirac 
particles. Numerical values of this correction are presented in the text. 
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The contribution to a.z{mz) from the splitting of 45h is found to be: 



Aa 3 (m z ) 4 5 H 



14tt 



[61n( 



M(l,l,±l) 



) - 241n( 



M(l,3,0) 
Mu 



+ 211n( 



M(8,l,0) 



)] (119) 



Here the numbers (a,b,c) in M denote the quantum numbers of the respective sub-multiplet 
in 45h under SU(3) x SU(2) x U{1). Within the simplest model consisting of a 45h only 
(i.e., if we ignore its couplings to other multiplets), one can argue that the masses of the sub- 
multiplets are about two (to one) order of magnitude lower than M\j. This is because with 
only a 45h, the effective SO (lO)-invariant superpotential including self-couplings has the 
form W 45lI = Mx(45h) 2 + k(45h) 4 /M where M characterizes the effect of quantum gravity. 
Note no cubic term (45 H ) 3 is allowed. Setting F 45h = M 1 < 45 H > +k < 45 H > 3 /M = 0, 
we get Mi ~ k < 45h > 2 /M « /.• \ / \ I . where we put < 45h >~ Mu. Putting k ~ 1 and 
M ~ Mpi = 2 x 10 18 GeV, we get M x ~ 10" 2 Mu. With the masses of the sub-multiplets 
given by M 1 w 10- 2 Mu, Eq. (119) yields Aa 3 (m z ) w -0.0045. 

Note that if one had introduced larger Higgs multiplets Ike 126h or even 54 of 5*0(10), 
one would have encountered substantially larger threshold correction to a 3 (m z )- For 
example, with just a 54, which contains a (6,1,4/3) + (6,1,-4/3) and a (8,1,0) sub- 
multiplets, the threshold correction from the color (sextets, octet) to a 3 (mz) is given by 
Aa 3 (m^) 6+ Q 8 = [(51, 21)/(147r)]a 3 (mz) 2 ln(M 6j8 /Mu). This yields a correction to a 3 {m z ) 
exceeding about 14% (of either sign) if M 6 8 /Mu ~ 1/2 or 2. It therefore appears that S'O(IO) 
multiplets with dimension larger than 45 are highly disfavored, if the observed unification 
of gauge couplings is not to be regarded as a coincidence. 
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